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1.  FORWARD 


Continuum  visible  emission  (the  visible  shuttle  glow)  is  observed  around  Space  Shuttle 
and  other  satellites  operating  in  low  earth  orbit.  The  emission  is  found  to  be  proportional  to  the 
local  0-atom  concentration  and  is  apparently  the  result  of  reaction  between  ambient  oxygen 
atoms  and  surface  adsorbed  NO.  The  detailed  kinetic  process  by  which  this  occurs  is  not  yet  well 
established.  The  SKIPPER  mission,  designed  to  monitor  the  bow  region  emissions  of  a  blunt 
vehicle  as  it  re-enters  from  low  earth  orbit,  provided  an  exceptional  opportunity  to  amass  critical 
data  on  this  phenomenon  which  could  have  been  used  to  validate  a  proposed  mechanism. 

The  potential  kinetic  processes  which  could  contribute  to  this  process  are  manifold.  They 
include  both  O  and  NO  surface  accommodation  and  chemical  reactivity,  collisional  and  thermal 
desorption,  and  various  non-Maxwellian  gas  phase  kinetic  interactions.  In  an  earlier  effort,  a 
detailed  glow  model  was  developed  and  the  available  literature  reviewed  to  provide  estimates  of 
the  various  required  kinetic  parameters.  Limited  experimental  evaluation  of  NO  accommodation 
and  collisional  desorption  coefficients  were  also  provided. 

The  present  research  effort  was  directed  towards  understanding  and  evaluating  the  kinetic 
mechanisms  responsible  for  the  visible  shuttle  glow.  The  technical  work  is  provided  in  three 
sections  as  described  below. 

•  An  evaluation  of  adsorption/desorption  behavior  of  NO  on  engineering  surfaces. 

It  is  believed  that  collisions  between  ambient  oxygen  atoms  and  surface  adsorbed  NO 
produce  the  visible  "shuttle  glow."  We  have  provided  measurements  which  can  be 
interpreted  in  terms  of  the  NO  surface  accommodation  coefficient  and  the  thermal 
desorption  coefficient.  These  are  key  parameters  in  the  kinetic  model  for  the  surface 
glow. 

•  A  critical  review  of  laboratory  measurements  of  visible  shuttle  glow. 

A  variety  of  papers,  some  seemingly  inconsistent,  have  been  published  on  laboratory 
studies  of  the  glow.  These  were  critically  reviewed  and  it  was  found  that  the  weight  of 
data  suggests  that  the. glow  occurs  through  a  combined  Eley-Rideal,  Langmuir- 
Hinshelwood  kinetic  mechanism.  An  interesting  data  collection  on  the  temperature 
dependence  of  the  glow  was  also  developed.  This  work  was  presented  as  an  invited  paper 
at  the  December  1994  AGU  meeting. 

•  A  collaboration  with  Iain  Boyd,  Cornell  University,  on  a  Monte  Carlo  model  for  the  glow. 

We  conclude  that  the  observed  altitude  dependence  of  the  glow  in  part  results  from 
transition  flow  effects.  A  DSMC  model  for  this  phenomena  has  been  assembled  and 
exercised. 

Several  reports,  publications  and  presentations  are  included  as  appendices. 
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2.  TECHNICAL  STUDIES 


2.1  NO  Adsorption/Desorption  Study 

The  surface  kinetics  responsible  for  the  shuttle  glow  can  be  quite  complex.  Knowledge  is 
required  on  the  surface  accommodation,  thermal  and  collisional  desorption  and  chemical  reaction 
of  chemical  species  such  as  N,  O  and  NO. 

A  few  years  ago,  we  reported  a  preliminary  study  of  NO  surface  accommodation  and 
thermal  desorption  (Caledonia  and  Krech,  1994).  There  was  some  uncertainty  in  those  earlier 
results,  and  in  the  present  program  we  have  repeated  and  expanded  the  study. 

Basically,  we  developed  a  simple  model  for  the  NO  surface  coverage  history  of  an 
initially  clean  material  in  vacuum  which  is  "instantaneously"  subjected  to  a  constant  background 
pressure  of  NO.  The  model  is  limited  to  monolayer  coverage;  the  experiment  was  performed  at 
room  temperature  where  the  NO  vapor  pressure  is  sufficiently  high  to  preclude  formation  of  a 
deeper  NO  surface  layer.  We  described  this  model  in  the  1994  report.  It  is  summarized  again 
below. 


Let  f  be  the  fractional  coverage  of  NO  on  the  surface,  S^q  be  the  sticking  coefficient,  and 
Tj,^0  be  the  characteristic  time  for  thermal  desorption.  Then 


^  =  (NO)  c 
dt  N  4 

S 


Sno  (1  -f) 


(2.1.1) 


where  N^  is  the  previously  defined  monolayer  surface  coverage,  (NO)  is  the  local  NO  concentra¬ 
tion  in  cm'^,  and  c  is  the  mean  NO  molecular  speed. 


This  has  the  steady  state  solution 

^  ^  S^o(NO)c/4 

S^o(NO)c/4  +  N, 


(2.1.2) 


which  can  be  recast  as 


5no(NO)  /- 


1-f  = 


N 


(0/4)1 


NO 


(2.1.3) 


Thus,  if  one  can  measure  f^  versus  NO  concentration,  one  could  measure  the  product  S^q  t^q. 
This  is  a  straightforward  steady-state  measurement. 
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The  time  dependent  solution  of  Eq.  (2.1.1)  provides  additional  information.  This  is 


/ 

f  =  fjl 


e 


(NO)  £ 
Ns  4 


%0 


(2.1.4) 


Thus,  by  monitoring  f  at  early  times,  after  the  NO  flow  is  turned  on,  one  can  back  out  the  overall 
characteristic  time,  whichever  dominates. 

In  our  experiment,  we  subjected  the  surface  to  a  controlled  NO  pressure  and  then 
bombarded  the  surface  with  8  km/s  oxygen  atoms  at  various  times  after  NO  exposure.  The  O- 
bombardment  results  in  NO2  fluorescence  where  the  fluorescence  intensity  is  necessarily 
proportional  to  the  surface  coverage  of  NO.  The  experimental  procedure  was  described  in  our 
1994  report  and  is  again  summarized  below. 

The  experimental  measurements  were  performed  within  our  0-atom  facility  on  a  target 
plate  15.2  cm  high  by  35.6  cm  long.  A  schematic  of  the  experiment  is  shown  in  Figure  2.1.1. 
The  plate  was  placed  some  72.5  cm  from  the  nozzle  throat  and  was  completely  engulfed  by  the 
0-beam  which  was  typically  pulsed  at  2  Hz.  A  typical  0-beam  fluence  per  pulse  at  the  plate  was 
3.6  X  lO^'^  cm'^,  a  fraction  of  a  monolayer.  The  plate  temperature  was  held  at  290  K. 


Source 

Chamber:  16  m.  O.D.  Tubing 

Six-way  Cross 
8  in.  O.D.  Tubing 

Puised  ^ 

Valve  and  [ 

Nozzle  -►DO 
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-  Target 

-  Laser 
Focussing 
Lens 


CO2 

Laser 


Gate 

Valve 


Cryogenic 

Pump 


B-1828a 


Figure  2.1.1  Schematic  of  experimental  apparatus. 
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The  glow  intensity  was  measured  just  above  the  plate  at  a  central  wavelength  of  630  nm 
with  a  10  nm  bandpass.  This  wavelength  is  just  to  the  blue  of  peak  glow  intensity.  The  optical 
field  of  view  was  cylindrical  with  a  3.2  cm  diameter. 

The  0-beam  exhibits  a  small  intrinsic  radiation  level  at  this  wavelength.  Furthermore  the 
beam  interacts  with  the  gaseous  NO  to  produce  a  continuum  glow.  This  gas  phase  glow  is  blue- 
shifted  related  to  the  surface  glow,  and  most  likely  results  from  energetic  oxygen  atom  inter¬ 
actions  with  dimers.  The  measurements  of  glow  intensity  above  the  surface  will  necessarily 
include  a  contribution  from  this  gas  phase  glow  as  well  as  from  the  surface  glow.  These  relative 
contributions  are  distinguished  by  performing  the  intensity  measurements  with  and  without  the 
target  in  place.  The  difference  in  intensities  between  these  two  measurements  is  taken  to  be  the 
contribution  of  the  surface  glow.  An  example  of  such  measurement  is  shown  in  Figure  2.1.2  for 
an  NO  pressure  of  3.7  x  10'^  torr.  As  can  be  seen,  the  gas  glow  intensity  is  small  compared  to 
that  of  the  surface  glow.  Note  that  the  initial  NO  exposure  occurred  at  slightly  different  times  for 
the  two  traces  shown  in  Figure  2.1.2. 


Figure  2. 1 .2  Pulse  rate  2  Hz,  measured  O  +  NO  surface  and  gas  glow.  P^q  =  3.7  x  10^ 

O  +  NO  surface  glow.  Bottom:  O  +  NO  gas  glow 

The  summary  of  "steady"  intensity  data  reported  in  1994  is  shown  in  Figure  2.1.3,  which 
is  taken  from  that  source.  The  data  shown  versus  pressure  is  the  saturated  surface  glow , 
concomitant  gas  phase  glow  and  the  difference  between  the  two  signals. 

As  pointed  out  earlier,  the  50%  saturation  point  in  the  data  occurs  at  ~10  ^  torr  and  thus 
fromEq.  (2.1.3) 


1 


=  3.3  X  10** 


^NO  ^NO 


4N, 


(2.1.5) 


4 


2.00 


1.80 

1.60 


1.40 

1  1.20 

o  1.00 
> 

c5 

(D  0.80 

cr 

0.60 


0.40 

0.20 


0.00 

10-6  10-5  10-4  10-3 

NO  Background  Pressure  (Torr) 

C-6494 

Figure  2.1.3  Observed  steady  glow  intensity  at  630  nm  and  290  K  versus  no  pressure,  8  km/s  O- 
atom  at  2  Hz  on  anodized  Al. 
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where  we  have  substituted  in  the  NO  concentration  at  P  =  10"^  torr.  for  c  =  4.5  x  10"^  and 
Nj  ~  3  X  10^^  cm‘^  we  have  then 


S 


NO’^'HO 


0.81s 


(2.1.6) 


and,  since  Sj,^q  cannot  exceed  unity,  this  implies  that  at  290  K 

Tno  ^  0.81s  . 

We  have  included  Eq.  (2. 1 .2),  evaluated  with  Expression  (2. 1 .6)  on  Figure  2.1.3.  As  can 
be  seen,  the  comparison  is  reasonable.  Nonetheless,  this  evaluation  is  uncertain  to  at  least  a  . 
factor  of  two  given  the  scatter  in  the  data  and  the  arbitrariness  in  the  specification  of  N^  (more 
properly,  Eq.  (2.1.3)  should  be  used  to  define  the  quantity  Sj..jq  %o^s)- 

We  have  now  extended  this  work  to  exploit  additional  information  which  can  be  found  in 
the  time  variation  of  the  glow  intensity.  More  specifically,  we  monitored  the  characteristic  e-fold 
time  of  the  intensity  at  NO  pressures  of  1.4  x  10'^,  2.3  x  10'^  and  3.7  x  10'^  torr.  These 
measurements  were  made  several  times  and  the  results  were  highly  repeatable.  Furthermore  at 
these  pressures,  the  gas  phase  glow  intensity  never  exceeded  20%  of  the  total  observed  intensity. 
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In  all  cases,  we  found  that  the  characteristic  e-fold  time  varied  inversely  with  pressure  as 
expected  from  Eq.  (2.1.5).  Indeed,  the  characteristic  e-fold  time  dropped  from  4.5  to  2  s  over  the 
aforementioned  pressure  range.  This  temporal  variation  is  slightly  less  than  linear  with  pressure 
presuming  reflecting  the  contribution  of  the  tj^jQterm  in  the  exponential  of  Eq.  (2. 1 .5).  We  can 
use  the  lowest  measured  e-fold  time  of  2  s  to  make  a  lower  bound  estimate  for  Sj^q 
same  values  mentioned  previously  for  c  and  N^.  From  this,  we  deduce  that  Sj^q  ^0.1  whence 
<  8.1  s.  These  values  are  of  course  consistent  with  the  50%  saturation  point  occurring  at  an 
NO  pressure  of  10'^  torr. 

Note  that  the  additional  measurements  of  the  e-fold  time  for  steady  coverage  allows  the 
uncertainty  in  specifying  N^  to  only  affect  the  uncertainty  in  Sj,jq,  i.e.,  a  factor  of  two  increase  in 
Ng  will  produce  a  factor  of  two  increase  in  Sj^q-  "^he  remaining  uncertainty  in  Tj^jq  is  only  that 
resulting  from  the  scatter  in  the  data. 

We  have  found  no  measurements  in  the  literature  for  comparison  with  the  present  results. 
2.2  A  Review  of  Laboratory  Studies  of  the  Visible  Shuttle  Glow 
2.2.1  Introduction 


There  have  been  many  studies  of  the  visible  Shuttle  glow  both  in  flight  (see  the  review  by 
Garrett  et  al.  (1988),  as  well  as  more  recent  measurements  by  Viereck  et  al.  (1991),  Viereck  et  al. 
(1992),  and  Swenson  et  al.  (1995)),  and  in  the  laboratory  (Arnold  and  Coleman,  1988,  1991; 
Caledonia  et  al.,  1990, 1993;  Orient  et  al.,  1990, 1992;  Swenson  et  al.,  1991;  Ardebili  et  al., 

1991;  and  Greer  et  al.,  1993).  It  is  now  well  established  that  the  visible  Shuttle  glow  arises  from 
the  surface  catalyzed  recombination  of  NO  and  O  forming  electronically  excited  NO2.  The 
mechanism  by  which  this  occurs  remains  elusive,  however.  In  this  work  we  will  review  the 
available  laboratory  data  on  surface  catalytic  recombination  of  NO2  and  identify  those  key 
observations  which  can  lead  to  a  consistent  interpretation  of  the  glow  mechanism. 

2.2.2  Review  of  Experimental  Data 

The  three-body  gas  phase  recombination  reaction  between  thermal  O  and  NO  has  been 
well  studied.  The  electronically  excited  NO2  produced  in  the  reaction  emits  a  yellow-green 
continuum  spectrum  known  as  the  air  afterglow.  A  representative  spectral  measurement 
(Fontijn  et  al.,  1964),  among  many,  is  shown  in  Figure  2.2.1.  When  similar  measurements  are 
performed  in  a  flow  of  O  and  NO  directed  over  a  nickel  surface,  an  orange  glow  is  observed,  see 
Figure  2.2.1  (Kenner  and  Ogryzlo,  1984;  Chu  et  al.,  1986;  and  references  therein). 

The  red  shift  of  the  spectrum  is  usually  presumed  to  be  due  to  energy  loss  to  the  surface 
in  the  heterogeneous  reaction.  Nonetheless,  although  this  reaction  has  been  studied  since  the 
1960's,  the  kinetic  interpretation  of  these  spectra  is  still  uncertain  due  to  the  possible  occurrence 
of  simultaneous  gas  and  surface  phase  processes.  Specifically,  these  measurements  are  made  in 
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1.0 


Figure  2.2. 1  Thermal  studies  of  O  +  NO  recombination  glow  contrasted  with  the  shuttle  glow, 

_ gas  phase  glow,  Fontijn  et  al.  (1964); _ shuttle  glow,  Viereck  et  al.  (1992); 

surface  catalyzed  glows; ....  Kenner  and  Ogryzlo  (1984); - ,  Chu  et  al.  (1986). 

mixed  streams  of  O  and  NO  at  pressures  exceeding  0.02  torr.  While  surface  catalyzed 
recombination  reactions  such  as: 


NO  +  O  A  NO2* 


(2.2.1) 


O  +  O  A  02(aX),  o; 


(2.2.2) 


will  occur,  the  pressures  are  sufficiently  high  such  that  subsequent  gas  phase  reactions  such  as: 

O3*  +  NO  NO2*  +  O2  (2.2.3) 

02(A)  +  O2  O3*  +  O  (2.2.4) 


can  also  contribute  to  the  observed  glows  (Kenner  and  Ogryzlo,  1984;  and  Chu  et  al.,  1986). 
Indeed,  Chu  et  al.  (1986)  demonstrated  that  the  effective  lifetime  of  the  glow,  as  determined  by 
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its  spatial  extent,  decreased  by  an  order  of  magnitude  when  O2  was  eliminated  from  the  flow, 
clearly  demonstrating  the  importance  of  gas  phase  reactions  subsequent  to  surface  recombination. 

Shown  for  comparison  with  the  results  of  these  thermal  studies  is  a  representative 
measurement  of  the  Shuttle  glow  (Viereck  et  al.,  1992).  This  glow  is  clearly  red-shifted  relative 
to  the  gas  phase  continuum  yet  is  to  the  blue  of  the  thermal  surface  catalytic  glows. 

There  have  now  been  several  measurements  of  the  visible  glow  performed  with  beams  of 
atomic  oxygen  at  pressures  low  enough  to  eliminate  the  possibility  of  significant  gas  phase  reac¬ 
tion.  Three  such  studies  (Arnold  and  Coleman,  1990, 1992;  Ardebili  et  al.,  1991;  Greer  et  al., 
1993)  have  been  performed  at  0-beam  energies  of  eV,  typical  of  arc-type  oxygen  atom 
sources.  In  these  works  a  fast  0-beam  (diluted  in  a  rare  gas),  flux  <3  x  10*^  cm  ^‘S  \  was  mixed 
with  an  effusive  source  of  NO  above  a  surface  (nickel,  metals,  Z306)  with  total  pressures  below 
10'^  torr. 

The  results  of  these  measurements,  contrasted  with  the  Shuttle  glow,  are  shown  in 
Figure  2.2.2.  Note  the  comparison  between  the  three  sets  of  measurement  is  quite  good 
(although  the  spectral  resolution  of  the  data  of  Greer  et  al.  (1993)  is  quite  coarse,  150  nm 
FWHM).  No  significant  variations  in  spectral  shape  were  observed  as  the  target  material  was 
changed. 


Figure  2.2.2  Low  velocity  O  beam  recombination  glows  contrasted  with  shuttle  glow. 

_ ,  Shuttle  glow,  Viereck  et  al.  (1992); ....  Ardebili  et  al.  (1991), 

Vq  =  1.85  km/s; -  Arnold  and  Coleman  (1990),  Vq  =  1.4  km/s;  symbols, 

Greer  et  al.  (1993),  ■  -Vq  =  1.5  km/s,  A  -  Vq  =  3.3  km/s. 
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There  are  at  least  three  significant  aspects  to  these  measurements: 

1 .  the  pressures  are  so  low  that  only  surface  reactions  should  pertain 

2.  the  measurements  of  Greer  et  al.  (1993)  show  no  0-atom  velocity  dependence  on  the 
glow  intensity 

3.  the  observed  spectrum  is  red-shifted  relative  to  the  thermal  atom  surface  catalyzed 
glows. 

This  last  point  is  particularly  interesting  in  that,  if  anything,  one  might  have  conjectured  a  blue 
shift  to  result  from  the  increased  O-atom  velocity.  Perhaps  this  suggests  that  a  different  kinetic 
mechanism  is  operative  in  the  free  molecular  regime.  We  will  return  to  this  point. 

Two  studies  of  the  surface  catalyzed  O  +  NO  recombination  glow  have  been  performed 
under  free  molecular  conditions  with  an  oxygen  beam  velocity  of  ~  8  km/s,  or  orbital  velocity.  In 
both  instances  the  beams  were  formed  from  pure  O2  without  a  rare  gas  diluent.  In  the  first  of 
these  (Caledonia  et  al.,  1990;  Swenson  et  al.,  1991),  a  pulsed  0-beam  was  used  to  irradiate  a 
surface  which  had  previously  been  dosed  with  less  than  a  mono-layer  of  NO.  The  0-beam  flux 
was  <  10^^  cm‘^«s'^  and  the  operating  pressure  <  10"^  torr.  In  the  second  experiment  (Orient  et 
al.,  1990,  1992),  a  CW  0-beam  with  flux  <  10^^  cm'^  s"^  and  an  NO  gas  jet  were  directed  onto 
surfaces  at  total  pressures  <2x10"^  torr.  In  both  experiments  targets  of  several  different 
materieils  were  used  with  no  significant  change  observed  in  the  spectral  shape  of  the  glow. 

The  results  of  these  two  measurements,  again  contrasted  with  the  Shuttle  glow,  are  shown 
in  Figure  2.2.3.  The  disparity  between  the  two  measurements  is  not  understood,  however  it  is 
clear  that  the  glow  induced  by  8  km/s  oxygen  atoms  is  significantly  blue  shifted  relative  to  that 
observed  with  lower  velocity  0-beams,  as  shown  in  Figure  2.2.2.  Furthermore,  the  measurement 
of  Caledonia  et  al.  (1990)  is  in  excellent  agreement  with  the  observed  Shuttle  glow  spectra, 
demonstrating  that  the  Shuttle  glow  can  be  simulated  in  the  laboratory  by  the  interaction  of 
8  km/s  oxygen  atoms  with  absorbed  NO.  It  was  further  demonstrated  (Swenson  et  al.,  1991)  that 
a  similar  simulation  relative  to  the  spatial  extent  of  the  glow  was  also  achieved. 

There  have  also  been, several  laboratory  observations  on  the  temperature  dependence  of 
the  glow.  The  low  velocity  beam  measurements  (Arnold  and  Coleman  1988,  1990;  Ardebili  et 
al.,  1991)  report  confusing  variations  of  glow  intensity  with  temperature.  These  include  one 
hundred-fold  changes  in  intensity,  apparently  due  to  surface  passivation  effects,  and  convoluted 
temperature  variations.  No  attempt  will  be  made  to  rationalize  these  observations  here.  On  the 
other  hand,  measurements  performed  with  the  8  km/s  beams  provide  an  interesting  and 
unexpected  comparison  with  each  other  and  with  flight  data  as  described  below. 

Swenson  et  al.  (1986)  examined  glow  intensities  from  several  shuttle  flights  and 
concluded  that  the  glow  intensity  exhibited  a  temperature  dependence.  They  proposed  an 
Arrhenius  type  dependence,  exp  AE/kT,  with  characteristic  energy  AE  -0.12  eV.  More  recently. 
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Figure  2.2.3 


8  km/s  0-beam  recombination  glows  contrasted  with  the  shuttle  glow, ....  Orient 
et  al.  (1992); - Caledonia  et  al.  (1990);  Shuttle  glow,  Viereck  et  al.  (1992). 


Orient  et  al.  (1992)  reported  laboratory  measurements  of  the  temperature  dependent  glow 
intensity  above  MgF2,  Ni  and  Ti  substrates  using  their  CW  8  km/s  O-beam  source.  Their  results 
which  cover  the  temperature  range  of  240  to  340  K  are  shown  in  Figure  2.2.4.  They  also  fit  their 
data  to  an  Arrhenius  expression,  deducing  values  of  AE  «0.1 1  to  0.12  eV  in  apparent  agreement 
with  the  suggestion  of  Swenson  et  al.  (1986). 

The  temperature  variation  of  the  glow  intensity  was  also  observed  with  a  pulsed  8  km/s 
O-beam  (Oakes  et  al.,  1992).  These  results  are  also  shown  in  Figure  2.2.4.  These  measurements 
were  performed  with  a  constant  NO  background  pressure  of  10'^  torr  while  the  target  temperature 
was  varied.  (The  data  may  be  suspect  below  150  K  because  of  condensation  on  the  substrate.) 
Clearly,  the  observed  variation  in  glow  intensity  with  temperature,  does  not  correlate  with  a 
simple  Arrhenius  behavior.  Nevertheless,  both  sets  of  data  are  in  remarkably  good  agreement 
over  the  common  temperature  range  of  the  measurements. 

This  comparison  suggests  that  an  Arrhenius  expression  is  not  appropriate  for  modeling 
the  glow  intensity.  Indeed,  as  can  be  seen  in  Figure  2.2.4,  an  Arrhenius  expression  would  require 
a  temperature  variation  which  would  strongly  diverge  from  the  measurement  of  Oakes  et  al. 
(1992)  at  lower  temperatures.  The  good  comparison  between  the  two  sets  of  data  suggest  that  the 
glow  intensity  temperature  dependence  is  more  complex  than  a  simple  activation  energy.  Most 
compelling,  very  recent  flight  data  reported  by  Swenson  et  al.  (1995)  arbitrarily  normalized  in 
Figure  2.2.4,  exhibit  a  weaker  temperature  dependence  than  predicted  by  the  earlier  Arrhenius 
expression  and  are  in  reasonable  agreement  with  the  laboratory  data. 
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Figure  2.2.4  Observed  temperature  dependence  of  O  +  NO  surface  glow  intensity,  5  eV 

0-atoms. _ Oakes  et  al.  (1992),  A  =  400  to  800  nm,  Al; ....  Orient  et  al. 

(1992),  X  =  625  nm,  MgF2,  Ni,  Ti;  contrasted  with  recent  shuttle  observation 

- Swenson  et  al.  (1995).  Also  shown  for  comparison  is  the  Arrhenius 

expression  exp  (0.12  eV/kT)  arbitrarily  normalized, - . 

A  last  point  with  regard  to  the  laboratory  measurements  concerns  the  production  of  NO. 
Various  mechanisms  have  been  proposed  for  NO  production  around  shuttle  (e.g.,  see  Swenson  et 
al.,  1985;  Green  et  al.,  1986;  Kofsky  and  Barrett,  1986).  These  include  surface  recombination  of 
ambient  N  and  O,  surface  dissociation  of  ram  N2  producing  N  which  can  ultimately  form  NO, 
ambient  NO,  and  the  gas  phase  atom  exchange  reaction: 

O  +  N2  NO  +  N  (2.2.5) 


which  is  thermoneutral  at  0-atom  ram  velocities.  Although  the  specific  NO  formation 
mechanism  remains  to  be  defined,  this  last  reaction  has  been  studied  in  the  laboratory  for  O-atom 
velocities  of  7  to  12  km/s  (Upschulte  et  al,  1992)  and  its  cross  section  has  been  shown  to  be 
large  enough  to  produce  NO  densities  consistent  with  observed  flight  NO2*  emission  levels.  The 
NO  formation  mechanism  is  particularly  important  in  that  it  helps  determine  the  scaling  of  glow 
intensity  with  atmospheric  species  concentration. 

Analysis 

We  are  in  the  process  of  developing  a  detailed  kinetic  mechanism  for  the  visible  glow 
which  will  allow  us  to  scale  the  glow  observations  to  other  flight  conditions.  In  this  section  we 
will  review  the  status  of  that  analysis  and  the  impact  of  the  laboratory  observations  on 
mechanistic  interpretation. 
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Surface  catalyzed  recombination  is  usually  considered  in  terms  of  two  limiting 
mechanisms.  The  first  of  these  is  the  Langmuir-Hinshelwood  (L-H)  where  recombination  is 
assumed  to  occur  between  two  adsorbed  species  with  subsequent  release  of  the  recombined 
species  into  the  gas  phase,  e.g. 


A(S)  +  B(S)  A  AB  * 


(2.2.6) 


where  the  symbol  S  refers  to  surface  bound  species.  Alternatively  the  Eley-Rideal  mechanism 
proposes  recombination  occurring  in  the  collision  between  an  incident  gaseous  species  and  an 
absorbed  species  with  ultimate  release  of  the  recombined  species  in  the  gas  phase,  e.g. 

A  +  B(S)  A  AB  *  .  (2.2.7) 


In  the  past,  either  one  or  the  other  of  these  mechanisms  has  been  proposed  to  interpret  the  visible 
shuttle  glow.  Ardebili  et  al.  (1991)  has  suggested  a  two  step  process  involving  both  mechanisms 
to  explain  the  complex  temperature  dependence  of  their  laboratory  observations. 

Assuming  that  the  source  of  NO  is  in  the  gas  phase,  e.g.  the  exchange  reaction  (2.2.5),  or 
even  ambient  NO,  then  a  number  of  kinetic  steps  could  be  involved  in  the  glow  process.  These 
are  summarized  in  Table  2.2.1 .  One  must  consider  surface  adsorption  and  collisional  and 
thermal  desorption  reactions  for  both  O  and  NO  as  well  as  the  ultimate  surface  recombination 
reactions.  If  the  Eley-Rideal  mechanism  pertains,  the  NO2  collisional  formation  mechanisms 
become  important  as  well.  Note  furthermore  that  the  NO2  production  reactions  will  not 
necessarily  produce  NO2  in  the  excited  state  with  unit  efficiency.  The  inclusion  of  all  these 
mechanisms  along  with  the  NO  production  mechanism  produces  a  necessarily  complex 
relationship  for  the  glow  intensity.  This  is  particularly  true  since  all  these  kinetic  quantities  can 
depend  on  both  temperature  and  material  (thus  the  importance  of  the  data  in  Figure  2.2.4).  Some 
limited  data  on  the  individual  mechanisms  is  available  for  engineering  material  surfaces.  For 
instance  collisional  desorption  efficiencies  of  adsorbed  NO  impacted  by  8  km/s  oxygen  atoms 
have  been  reported  (Krech  at  al.,  1993)  as  0.05  to  0.01  for  selected  materials.  As  noted  in 
Section  2.1,  we  have  found  room  temperature  NO  sticking  coefficients  ^0.1  with  thermal 
desorption  times  of  less  than  8  seconds.  Ultimately,  the  available  visible  laboratory  and  flight 
data  provide  the  critical  test  for  mechanistic  assumptions/interpretations. 

As  a  specific  example  let  us  review  the  implications  of  the  wavelength  dependence  of  the 
glow  observations  reported  in  the  previous  section.  The  glow  from  lower  energy  (^1  eV)  beam 
measurements  (Arnold  and  Coleman  (1990;  1992),  Ardebili  et  al  (1991),  Greer  et  al.  (1993))  are 
shifted  to  the  red  of  those  from  thermal  measurements  performed  at  higher  pressures  even 
though,  in  at  least  one  instance,  the  free  molecular  directed  velocity  was  close  to  thermal.  One 
interpretation  of  this  observation  is  that  the  free  molecular  measurements  represent  the  true 
surface  reaction  between  O  and  NO  while  the  higher  pressure  results  reflect  contributions  from 
other  species  and/or  gas  phase  reactions. 
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Table  2.2. 1  Suggested  Kinetic  Processes  of  Relevance  to  the  Visible  Shuttle  Glow  Mechanism 


Reaction 

Type 

Process 

Kinetic 

Quantity 

Adsorption 

0(g) - >0(S) 

NO(g) - »NO(S) 

Slicking 

coefficient 

Thermal 

Desorption 

0(S) - >0(g) 

NO(S) - >NO(g) 

Thermal 
desorption  rate 

Collisional 
desorption 
or  reaction 

0(g)+0(S) - »20(g) 

- >02(9) 

X(g)+0(S) - >X(g)+0(g) 

0(g)+N0(S) - yO(g)+NO(g) 

- »N02(g) 

0(g)+N02(S)^0(g)+N02(g) 

N0(g)+0(S)-^N0(g)+0(g) 

^N02(g) 

X(g)+NO(S) - »X(g)+NO(g) 

Collisional 
desorption 
or  reaction 
efficiency 

Surface 

Reaction 

0(S)  +  N0(S) - >N02(g). 

N02(S) 

0(S)+0(S) - >02(g).02(S) 

Reaction  rate 

*g  =>  Gas  Phase  Species 
*S  ==>  Surace  Adsorbed  Species 


On  the  other  hand  the  blue  shift  observed  in  the  higher  energy  (~5  eV)  beam 
measurements  (Orient  et  al.  (1990;  1992),  Caledonia  et  al.  (1990),  Swenson  et  al.  (1991))  is 
difficult  to  reconciliate  with  an  L-H  mechanism.  It  is  certainly  not  consistent  with  complete 
surface  accommodation  of  the  oxygen  atoms  prior  to  recombination,  although  presumably  a 
quasi-L-H  type  recombination  process  could  occur  with  partial  accommodation. 

Furthermore  Oakes  et  al.  (1992)  have  observed  a  velocity  squared  dependence  of  the 
glow  intensity  on  0-atom  velocity  over  the  velocity  range  of  6  to  12  km/s.  Again  this 
observation  is  more  consistent  with  an  Eley-Rideal  mechanism  than  an  L-H  mechanism.  Recall 
however,  that  Greer  et  al.  (1993)  did  not  observe  a  velocity  dependence  in  glow  intensity  at  the 
lower  O-atom  velocities  of  1.5  to  3.3  km/s. 

These  observations  collectively  suggest  that  the  visible  shuttle  glow  results  either  from  a 
Eley-Rideal  mechanism  or  more  likely  a  combination  mechanism.  One  possibility  would  be  L-H 
recombination  to  produce  free  and  bound  NO2,  e.g. 

0(S)  +  NO(S)  N02(g),  N02(S)  (2.2.8) 
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both  with  some  efficiency  for  producing  electronically  excited  NO2,  followed  by  a  knockoff  or 
exchange  reaction  with  incoming  oxygen  atoms,  e.g. 

0(g)  +  NO^CS)  0(g)+N02(g)  (2.2.9) 

again  with  some  efficiency  for  releasing  electronically  excited  gaseous  NO2. 

A  similar  mechanism  was  previously  proposed  by  Ardebili  et  al.  (1991)  to  rationalize 
their  experimental  observations. 

One  additional  measurement  which  supports  this  latter  hypothesis  has  been  reported  by 
Caledonia  et  al.  (1993).  In  this  work  surfaces  in  vacuum  were  dosed  with  NO2  and  then  irradi¬ 
ated  with  8  km/s  oxygen  atoms.  The  resulting  glow  was  spectrally  similar  to  that  observed 
when  surfaces  were  dosed  with  NO.  Although  this  observation  is  consistent  with  reaction  (2.2.9) 
it  must  be  cautioned  that  potential  effects  due  to  dissociative  chemisorption  of  NO2  cannot  be 
ruled  out. 

Summary 

In  summary,  the  data  base  on  laboratory  studies  of  surface  catalyzed  recombination  of 
O  +  NO  has  been  reviewed  and  analyzed  in  terms  of  the  visible  shuttle  glow.  Observed  trends 
in  the  laboratory  data  have  led  to  the  suggestion  of  a  combined  Eley-Rideal,  Langmuir- 
Hinshelwood  kinetic  mechanism  for  production  of  the  laboratory  glow.  The  scaling  behavior  of 
this  mechanism  is  presently  being  examined  in  the  light  of  flight  observations. 

It  should  be  noted  (Gelb,  1995)  that  all  laboratory  studies  to  date  have  involved  energetic 
oxygen  atoms  impacting  adsorbed  NO.  In  flight  the  inverse  reaction,  energetic  NO  impacting 
adsorbed  O,  can  also  occur  and  this  could  result  in  different  glow  observations. 

2.3  Altitude  Dependence  of  the  Glow  and  DSMC  Modeling 

There  have  only  been  two  measurements  of  the  altitude  variation  in  the  shuttle  glow 
intensities  taken  during  an  individual  flight  where  all  other  experimental  parameters  are  held 
fixed.  The  first  of  these  is  from  the  Atmosphere  Explorer  series  of  the  1970's  and  was  reported 
by  Yee  and  Abreu  in  1983.  Their  results,  at  a  wavelength  of  732  nm,  are  shown  in  Figure  2.3.1. 
The  glow  data  span  the  altitude  range  of  140  to  270  km  and  as  can  be  seen,  the  intensity  scales 
linearly  with  the  oxygen  atom  density  (simultaneously  measured)  from  270  to  160  km  and  then 
begins  to  increase  more  rapidly  with  decreasing  altitude,  approximately  trailing  the  N2  density 
variation. 


Figure  2.3.1  Visible  glow  intensity  at  7320  A  as  observed  on  Atmospheric  Explorer  contrasted 
with  measured  oxygen  atom  and  nitrogen  molecule  concentrations  (Yee  and 
Abreu,  1983). 

For  many  years,  this  was  the  only  altitude  dependent  data  available  and  the  change  in 
slope  had  been  largely  interpreted  in  terms  of  a  change  in  the  dominating  chemical  kinetic 
mechanism.  Very  recent  measurements  reported  by  Swenson  et  al.  (1996)  and  shown  in 
Figure  2.3.2  may  change  this  interpretation.  This  data,  all  taken  in  one  shuttle  flight,  STS-62, 
shows  glow  intensity  tracking  the  0-atom  density  down  to  an  altitude  of  ~220  km  and  then 
increasing  less  steeply  than  the  0-atom  density  at  lower  altitude,  i.e.,  opposite  from  the 
Atmospheric  Explorer  observation. 

Swenson  et  al.  (1996)  suggest  this  falloff  may  be  the  result  of  density  buildup  on  the  ram 
side  of  shuttle  where  the  measurement  was  made.  Interestingly,  the  AE  satellite  is  approximately 
an  order  of  magnitude  smaller  than  shuttle.  Since  the  atmospheric  scale  height  is  =  60  km  at 
these  altitudes,  the  BCnudsen  number  for  the  shuttle  at  220  km  is  the  same  as  that  for  Atmospheric 
Explorer  at  140  km,  i.e.,  the  observed  change  in  slope  of  glow  intensity  occurs  at  the  same 
Knudsen  number  (or  for  that  matter,  Reynolds  number)  on  both  flights.  This  observation  of 
aerodynamic  scaling  strongly  suggests  that  the  glow  slope  "kink"  results  from  transitional  flow 
rather  than  chemical  kinetic  effects. 
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Signal  Intensity 


C-8190 


Figure  2.3.2  Observed  variation  of  glow  intensity  on  the  Experimental  Investigation  of  Shuttle 
Glow  (EISG)  experiment  on  space  shuttle  (Swenson  et  al.,  1996). 


The  Knudsen  number,  which  is  the  ratio  of  the  mean  free  path  to  characteristic  body 
dimension,  is  ~25  at  the  glow  intensity  "kink"  altitude.  Although  this  magnitude  Knudsen 
number  is  generally  considered  to  be  well  into  the  free  molecular  limit,  Probstein  (1961)  pointed 
out  that  for  hypersonic  flow,  it  is  more  appropriate  to  apply  the  hyperthermal  mean  free  path 

K  "  (2.3.1) 

where  is  the  flow  Mach  number,  =25  for  our  interest.  The  hyperthermal  mean  free  path  so 
defined  reflects  the  collision  mean  free  path  between  incident  and  surface  reflected  molecules. 

Indeed,  Bird  (1990)  recently  studied  space  shuttle  data  on  lift  over  drag,  L/D,  with  his 
Direct  Simulation  Monte  Carlo  (DSMC)  model  and  deduced  that  the  shuttle  data  was  most 
consistent  with  diffuse  reflection  of  ambient  molecules  off  of  shuttle  surfaces.  A  comparison  of 
shuttle  data  to  Bird’s  DSMC  calculations  is  shown  in  Figure  2.3.3  contrasted  with  the  free 
molecular  limit.  Bird  states  that  transition  regime  effects  on  space  shuttle  persist  to  about 
300  km.  More  recently,  Dogra  and  Moss  (1991)  also  predict  significant  transition  flow  effects 
on  space  shuttle  above  200  km. 

Thus,  there  is  strong  evidence  that  the  aforementioned  glow  observations  are  impacted  by 
transitional  flow  effects  and  that  a  simple  free  molecular  model  is  insufficient  for  flight 
predictions.  The  fact  that  the  two  flights  show  differing  trends  for  glow  intensity  variation  at  the 
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Figure  2.3.3  Comparison  of  shuttle  measurements  of  lift  over  drag  with  DSMC  calculations 
assuming  diffuse  reflection  (Bird,  1990). 

lower  altitudes  is  intriguing.  Recall  that  the  glow  results  from  a  surface-catalyzed  reaction. 
Surface  properties  of  course  will  not  scale  with  Knudsen  number  and  will  be  altitude  dependent. 
The  surface  coverages  at  the  glow  "kink"  points  of  the  two  flights  should  be  significantly 
different  and  this  may  provide  for  the  difference  in  the  two  observations.  Again,  this  may  be 
addressed  in  a  more  sophisticated  model. 

We  have  entered  into  a  collaborative  effort  with  Prof.  Iain  Boyd,  Cornell  University,  to 
develop  transition  flow  predictions  of  visible  glow  intensity.  The  object  is  to  have  Prof.  Boyd 
exercise  his  DSMC  model  to  provide  surface  fluxes  and  velocity  distributions  of  the  key  species 
O,  N,  NO  and  N2  These  four  species  are  all  present  in  the  ambient  atmosphere,  but  N  and  NO 
can  also  be  produced  by  chemical  reaction  between  incident  and  scattered  O  and  N2,  as  described 
below.  These  flux  predictions  will  be  vehicle  specific,  i.e.,  scale  with  Knudsen  number,  and  can 
be  used  with  a  separate  surface  model  employing  relevant  kinetic  mechanisms  to  predict  the 
altitude  dependent  glow  intensity  behavior. 

The  key  to  the  DSMC  model  is  the  proper  calculation  of  the  impinging  flux  incorporating 
effects  due  to  surface  scattered  primary  ambient  species  O  and  N2  undergoing  secondary 
collisions  with  incoming  species,  including  chemical  reaction.  We  assume  the  key  chemical 
reaction  to  be  between  impinging  and  scattered  O  and  N2,  i.e., 

6  +  N2  -  NO  +  N  (2.3.2a) 

O  +  N2-NO  +  N  (2.3.2b) 
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Note  that  this  reaction  also  forms  N  atoms  which  can  also  reach  the  surface  and  ultimately  form 
NO2*. 


Reaction  2.3.2  is  thermoneutral  at  an  0-atom  velocity  of  ~8  km/s,  orbital  velocity,  so  for 
the  reaction  to  have  an  appreciable  rate  constant,  the  scattered  species  must  have  a  finite  velocity. 
Specifically,  the  center  of  mass  energy  of  the  collision  pair  must  exceed  the  reaction  endother- 
micity  of  3.26  eV.  The  reduced  mass,  M,  of  the  collision  pair  is  10.2  gms/mole.  The  collision 
pair  energy  in  the  center  of  mass  system  (which  automatically  conserves  momentum)  is 


where  Vj.gi  is  the  magnitude  of  the  relative  velocity  between  the  two  particles,  Vj  +  V2.  At 
8  km/s  relative  velocity  is  3.38  eV  for  the  0,N2  collision  pair.  The  cross  section  for  reaction 
increases  dramatically  with  Vj.gj.  (We  have  an  approximate  measurement  as  described  below.) 
Thus,  the  velocity  at  which  the  0,N2  scatters  from  the  surface  is  critical  to  predict  the  NO 
production  rate.  Unfortunately,  there  is  very  little  data  on  this  subject.  We  have  some  very 
preliminary  data  on  scattering  of  8  km/s  0-atoms  from  engineering  surfaces.  We  find  the 
scattering  is  isotropic,  but  that  the  scattered  velocity  is  not  thermal  but  of  order  1  km/s  with  a 
broad  spread.  Limited  data  on  N2  suggest  a  scattered  velocity  <1.6  km/s  at  low  energies.  This  is 
all  for  normal  incidence.  This  subject  is  a  research  area  in  itself.  We  have  suggested  that  for  our 
present  calculations  we  assume  isotropic  scattering  and  scattered  velocities  of  1  and  2  km/s. 
These  values  of  scattered  velocity  are  sufficiently  slow  that  the  center  of  mass  velocity  of  the 
collision  pair  is  towards  the  scattered  surface  (recall  this  discussion  is  limited  to  ram  flow). 

For  simplicity,  we  assume  the  NO,N  created  in  the  process  moves  towards  the  surface  at 
the  center  of  mass  velocity  (and  direction)  of  the  collision  pair 


M,  +M, 


when  i  and  s  signify  incident  and  scattered  and  once  again  1  and  2  can  signify  either  O  or  N2. 

We  will  neglect  an  angular  spread  on  the  reaction  products  for  the  present.  This  can  be  included 
at  a  later  date. 

There  are  no  data  on  the  cross  section  for  Reaction  (2.3.2).  Upschulte  et  al.  (1992)  report 
results  for  production  of  vibrationally  excited  NO  by  this  reaction.  These  measurements  are  for  a 
-1000  K  0-atom  collision  partner.  We’ve  extrapolated  total  cross  sections  for  Reaction  (2.3.2) 
from  this  data  to  be  2, 4,  6  and  15  x  10'^^  cm^  for  relative  collision  velocities  of  8,  9,  10  and 
1 1  km/s  respectively. 
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The  remaining  quantities  needed  are  spacecraft  dimensions.  The  AE  satellite  diameter  is 
1.35  m.  The  shuttle  is  more  problematic.  The  experiment  of  Swenson  et  al.  (1996)  involved 
measurement  across  a  1  x  1  m  structure.  Nonetheless,  the  stmcture  itself  was  in  the  bay  and  that 
general  environment  is  defined  by  the  orbiter  itself.  The  shuttle  body  is  ~  5  x  30  m  with  a  rear 
wing  spread  of  approximately  20  m. 

The  DSMC  modeling  is  in  process.  Preliminary  results  were  presented  at  the  June  1996 
AIAA  Thermophysics  Conference.  See  Appendix  A. 

Since  this  modeling  was  not  completed  in  time  for  the  Skipper  launch,  pre-flight 
predictions  were  provided  by  extrapolating  the  AE  data.  These  results  were  presented  in  Boyd  et 
al.  (1995)  and  are  included  as  Figure  2.3.4. 


Figure  2.3.4  Skipper  glow  intensity  predictions  at  732  nm. 
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A.  Reports 

1 .  G.E.  Caledonia  and  R.H.  Krech,  Technical  Progress  Report  for  1  June  through 
31  December  1994,  January  27,  1995. 

2.  ID.  Boyd,  D.P.  Karipides,  D.A.  Levin,  and  G.E.  Caledonia,  "Pre-Flight  Predictions 
of  Ultraviolet  and  Visible  Emissions  for  the  Skipper  Flight  Experiment,"  Cornell 
Univ.  Report  FDA-95-09,  December  1995. 

3.  D.P.  Karipides,  I.D.  Boyd,  and  G.E.  Caledonia,  "Development  of  a  Monte  Carlo 
Overlay  Method  with  Application  to  Spacecraft  Glow,"  AIAA-96-1847,  31st  AIAA 
Thermophysics  Conf.,  June  17-20,  1996,  New  Orleans,  LA. 

4.  G.E.  Caledonia,  R.H.  Krech  and  A.  Leone,  "A  Review  of  Laboratory  Studies  of 
Spacecraft  Atmospheric  Interaction  Glows,"  invited  paper  at  the  Fall  American 
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Abstract 

Glow  brightness  calculations  for  the  Atmospheric 
Explorer  satellite  are  made.  A  simple  model  for  glow 
is  proposed,  identifying  nitric  oxide  as  the  important 
species  in  the  gas  phase  controlling  glow  brightness. 
The  flow  field  is  determined  using  the  direct  Simula^ 
tion  Monte  Carlo  method.  A  novel  overlay  technique 
is  used  to  capture  the  microscopic  behavior  of  rare 
atmospheric  species,  in  particular  nitric  oxide.  Sen¬ 
sitivity  to  altitude  and  to  different  chemical  cross- 
sections  for  nitric  oxide  production  is  assessed.  The 
importance  of  nitric  oxide  production  is  found  to  be 
greatest  at  lower  altitudes.  At  higher  altitudes,  the 
ambient  concentration  of  nitric  oxide  is  the  critical 
factor  in  determining  glow  brightness. 

Introduction 

Ram  surfaces  in  low  earth  orbit  are  known  to  ex¬ 
hibit  visible  and  ultraviolet  glow.  Such  glows  have 
been  detected  during  space  shuttle  missions^ and 
satellite  flights  of  the  Atmospheric  Explorer  (AE).^ 
In  proposed  models  of  this  phenomenon,  species 
produced  through  chemistry  are  assumed  to  travel 
through  the  flow  field  without  undergoing  collisions 
and  gas-surface  reactions  are  entirely  responsible  for 
glow  production.^"^  These  models  neglect  two  im¬ 
portant  effects  that  are  critical  in  the  production  of 
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glow.  At  the  lower,  more  dense  altitudes,  a  weak 
shock  forms  in  front  of  the  ram  surface.  In  this 
region,  the  assumption  of  a  collisionless  flow  is  in¬ 
appropriate.  At  the  higher  altitudes,  production 
through  chemical  reactions  becomes  negligible  and 
the  inclusion  of  free  stream  concentrations  is  criti¬ 
cal.  From  data  taken  by  the  AE  satellite,  it  is  evi¬ 
dent  that  the  rate  at  which  the  brightness  at  6563  A 
changes  with  altitude,  and  therefore  with  free  stream 
density,  is  dramatically  different  below  160  km  com¬ 
pared  to  that  above  160  km.^  This  experimental  ev¬ 
idence  suggests  that  gas-gas  collisions  become  im¬ 
portant  below  160  km,  indicating  that  the  flow  is 
collisioncd.  Therefore,  to  predict  glow  brightness  ac¬ 
curately,  a  method  that  includes  the  effects  of  gas- 
gas  collisions  in  the  weak  shock,  ambient  species  con¬ 
centrations  and  chemical  production  throughout  the 
flow  field  is  required. 

The  free  stream  conditions  at  altitudes  around 
160  km  are  extremely  rarefied  and  a  high  degree  of 
both  thermal  and  chemical  non-equilibrium  is  to  be 
expected  throughout  the  flow  field.  The  direct  sim¬ 
ulation  Monte  Carlo  (DSMC)  technique  is  ideal  for 
simulating  rarefied  flows  and  has  been  successfully 
applied  to  hypersonic,  non-equilibrium  aerothermo- 
chemistry  in  previous  work.  Production  of  the  chem¬ 
ical  species  critical  to  the  formation  of  glow  in  this 
regime  is  controlled  by  rare  collisions  and  rare  chem¬ 
ical  events.  Proper  modeling  of  such  rare  events  us¬ 
ing  the  DSMC  technique  requires  modifications  to 
the  standard  technique.  These  changes  are  outlined 
in  this  paper.  In  particular,  procedures  for  handling 
the  interaction  of  rare  particles  with  a  more  dense 
background  flow  field  and  for  calculating  the  pro¬ 
duction  of  particles  due  to  rare  chemical  events  are 
described.  Results  are  presented  for  simulation  of 
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the  Atmospheric  Explorer  satellite  over  the  range  of 
altitudes  from  140  km  to  180  km.  General  flow  field 
structure,  various  velocity  distributions  and  calcula¬ 
tions  for  glow  brightness  are  discussed. 

Glow  Production  Model 

A  simple  model  for  the  glow  mechanism  is  proposed. 
The  production  of  glow  is  assumed  to  be  entirely  due 
to  the  gas-surface  reaction: 

N0{g)  +  0(a)  ^  NO;(g) 

NO2  NO2  +  hy 

where  (g)  indicates  a  molecule  in  the  gas  state 
and  (a)  indicates  a  surface  adsorbed  particle.  The 
NO2  is  formed  in  an  excited  state  and  then  sponta¬ 
neously  decays  to  the  ground  state,  emitting  a  pho¬ 
ton.  These  photons  comprise  the  glow.  The  molecu¬ 
lar  surface  coverage  of  a  satellite  results  from  a  com¬ 
plex  balance  of  atmospheric  and  flow  properties,  sur¬ 
face  accommodation,  thermal  desorption,  collisional 
desorption  and  chemical  reactions.  In  the  present 
study  the  simplifying  approximation  is  made  that 
the  surface  is  saturated  with  oxygen  atoms,  one  of 
the  dominant  ambient  species  at  the  altitudes  of  in¬ 
terest.  Obviously,  the  flux  of  nitric  oxide  striking  the 
surface  will  be  crucial  to  determining  the  brightness 
of  the  glow.  Where  the  flow  is  collisional,  gas-gas 
collisions  and  reactions  will  be  important  in  deter¬ 
mining  the  nitric  oxide  flux. 

A  number  of  simplifying  assumptions  are  made  in 
implementing  this  model.  Without  data  to  suggest 
otherwise,  every  NO  molecule  that  strikes  the  sur¬ 
face  is  assumed  to  interact  with  an  adsorbed  oxygen 
atom  and  the  probability  of  a  reaction  forming  NO2 
is  set  to  unity.  Clearly,  the  probability  of  reaction 
should  be  less  then  one.  Additional  modeling  is  pos¬ 
sible  that  would  account  for  more  detailed  behav¬ 
ior  of  the  surface  chemistry  and  production  of  NO2. 
Examples  of  such  potential  additions  to  the  model 
include  a  dependence  of  the  reaction  probability  on 
the  incident  energy  of  the  NO  molecule  and  a  vari¬ 
ation  of  the  surface  coverage  of  atomic  oxygen  with 
altitude.  Given  the  microscopic  detail  available  in  a 
DSMC  simulation,  the  framework  exists  to  readily 
incorporate  these  affects. 

Atmospheric  Considerations 

For  the  lower  range  of  altitudes  where  glow  bright¬ 
ness  data  was  taken  on  the  Atmospheric  Explorer 
(less  than  180  km),  the  ambient  concentrations  of 


NO  and  N  are  less  than  the  concentrations  of  N2, 
O2  and  O  by  several  orders  of  magnitude.  Further¬ 
more,  the  rates  of  the  relevant  chemical  reactions 
are  such  that  NO  and  N  are  not  produced  in  quan¬ 
tity  throughout  the  flow  field.  Using  the  definition 
that  any  species  with  a  mole  fraction  less  than  10~^ 
is  considered  rare  and  is  otherwise  considered  com¬ 
mon,  then  in  all  the  cases  considered,  NO  and  N  can 
be  considered  rare  species,  while  N2,  O2  and  O  are 
common  species.  Given  that  the  ratio  of  the  num¬ 
ber  density  of  the  common  species  to  that  of  the  rare 
species  is  so  great,  it  is  a  reasonable  assumption  that 
the  collisions  and  reactions  of  the  rare  species  do  not 
significantly  affect  the  common  species  flow  field. 

Flow  of  air  against  a  ram  surface  at  orbital  speeds 
(8  km/s)  is  primarily  a  five  species  gas  mixture, 
consisting  of  N2,  O2,  NO,  N  and  0.  Atmospheric 
concentrations  for  all  of  these  species  excluding  NO 
are  determined  using  the  MSIS-90  model. A  ro¬ 
bust  atmospheric  model,  dependent  on  a  wide  vari¬ 
ety  of  specific  input  conditions,  the  MSIS  model  is 
used  here  to  generate  a  generic  atmosphere.  Unfor¬ 
tunately,  the  number  density  of  nitric  oxide  is  not 
included  in  the  atmospheric  structure  produced  by 
the  model.  Determining  accurate  NO  number  den¬ 
sities  is  a  difiicult  problem.  The  formation  and  de¬ 
struction  of  nitric  oxide  is  a  key  component  of  the 
odd-nitrogen  cycle  that  occurs  in  the  Earth’s  atmo¬ 
sphere.  Many  different  factors  determine  the  con¬ 
centration  of  NO  at  any  particular  time  and  the 
concentration  is  highly  sensitive  to  many  of  these 
factors.  In  this  study,  estimates  for  nitric  oxide 
number  densities  are  taken  from  experimental  rocket 
measurements.^^  When  analyzing  the  simulation  re¬ 
sults,  it  should  be  kept  in  mind  that  only  estimates 
of  the  NO  concentrations  have  been  used.  These 
estimates  are  likely  different  from  the  actual  atmo¬ 
spheric  NO  concentrations  experienced  during  the 
AE  flights. 

DSMC  Resolution  Difficulties 

Simulating  rare  species  and  low  probability  events 
using  the  DSMC  technique  causes  a  number  of  reso¬ 
lution  difficulties.  When  simulating  a  flow  field  with 
one  or  more  rare  species,  there  is  often  no  accept¬ 
able  choice  for  the  value  of  the  particle  weight,  Wp, 
defined  as  the  ratio  of  the  number  of  real  particles 
to  the  number  of  simulated  particles.  A  relatively 
low  value  of  Wp  is  required  in  order  for  a  statis¬ 
tically  significant  number  of  particles  representing 
the  rare  species  to  be  generated.  However,  with  the 
particle  weight  set  at  such  a  value,  inordinate  num- 
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bers  of  particles  representing  the  common  species 
will  also  be  generated.  This  results  in  unacceptably 
high  numbers  of  simulated  particles  in  each  cell,  re¬ 
quiring  massive  amounts  of  computational  effort  to 
calculate  the  flow  field.  Alternatively,  the  value  of 
Wp  can  be  set  relatively  high,  so  that  reasonable 
numbers  of  common  species  particles  are  produced. 
In  this  case,  however,  very  few,  if  any,  rare  species 
particles  are  then  produced  and  the  simulation  no 
longer  generates  meaningful  statistics  for  the  rare 
species.  Prior  experience  with  rare  species  and  res¬ 
olution  difficulties  are  discussed  in  Refs.  12  and  13. 

Previous  attempts  to  handle  these  resolution 
difficulties  have  used  a  continuum-based  overlay 
method. A  two  step  procedure  is  used  to  calcu¬ 
late  the  flow  field.  In  the  first  step,  a  DSMC  simula¬ 
tion  is  performed  to  determine  the  properties  of  the 
common  species  throughout  the  flow  domain.  The 
particle  weight  in  this  simulation  is  set  so  that  es¬ 
sentially  no  rare  species  particles  are  created.  The 
second  step  involves  a  continuum-based  calculation, 
where  mass  conservation  equations  are  solved  for  the 
rare  species  using  the  underlying  DSMC  flow  field  so¬ 
lution.  The  assumption  is  made  that  the  behavior  of 
the  rare  species  does  not  have  an  appreciable  effect 
on  the  common  species  flow  field.  This  continuum- 
based  approach  suffers  from  two  problems.  The  pro¬ 
cedure  used  to  solve  the  mass  conservation  equations 
makes  use  of  the  temperature  calculated  for  the  com¬ 
mon  species  in  each  computational  cell.  However,  in 
the  range  of  altitudes  relevant  for  spacecraft  glow, 
the  flow  field  is  sufficiently  rarefied  and  in  such  a 
degree  of  non-equilibrium  that  thermodynamic  vari¬ 
ables  such  as  temperature  are  not  valid  representa¬ 
tions  of  the  state  of  the  gas.  The  velocity  distribu¬ 
tion  of  the  molecules  is  very  often  bimodal  and  cer¬ 
tainly  not  Maxwellian.  It  is  therefore  incorrect  to 
use  temperature  to  characterize  the  underlying  flow 
field.  In  addition  to  this,  the  use  of  a  continuum- 
based  overlay  also  results  in  a  loss  of  microscopic 
detail  of  the  rare  species.  Such  detail,  particularly 
velocity  distributions,  may  have  a  significant  impact 
on  chemical  rates  and  glow  production.  It  is  impor¬ 
tant  that  this  microscopic  information  is  retained 
throughout  the  simulation. 

DSMC  Overlay  Method 

In  this  study,  a  DSMC-based  overlay  technique  is 
developed.  Conceptually,  this  is  very  similar  to 
the  continuum  overlay  approach,  where  two  sepa¬ 
rate  simulations  are  performed  except  that  in  this 
case  both  simulations  use  the  DSMC  technique.  In 


the  first  {base)  simulation,  the  particle  weight  is  set 
so  that  the  common  species  are  represented  at  an 
appropriate  level.  In  the  second  (overlay)  simula¬ 
tion,  the  particle  weight  is  reduced  so  that  the  rare 
species  are  represented.  Generation  of  the  common 
species  is  suppressed  during  the  overlay  simulation. 
In  essence,  the  overlay  simulation  can  be  considered 
to  be  a  magnified  simulation,  focusing  only  on  the 
rare  species.  As  a  particle  method,  DSMC  is  ideal 
for  handling  the  non-equilibrium  nature  of  the  flow 
field  and  for  determining  the  microscopic  behavior 
of  the  particles.  Thus,  the  DSMC  overlay  approach 
overcomes  the  two  main  disadvantages  of  the  con¬ 
tinuum  overlay  method. 

Coupling  Between  Rare  and  Common  Species 

The  main  challenge  of  performing  a  DSMC  overlay 
simulation  lies  in  achieving  a  detailed  coupling  be¬ 
tween  the  two  distinct  simulations.  This  is  not  a 
concern  in  the  base  simulation  as  it  is  assumed  that 
the  flow  field  of  the  common  particles  is  unaffected 
by  the  behavior  of  the  rare  particles.  It  is  a  signifi¬ 
cant  concern,  however,  when  performing  the  overlay 
simulation.  The  overlay  simulation  cannot  simply  be 
assumed  to  be  a  regular  DSMC  simulation  at  a  lower 
density.  The  rare  particles  are  moving  through  and 
colliding  with  a  background  of  much  denser  common 
species.  In  addition,  rare  chemical  events  in  the  base 
simulation  that  did  not  noticeably  affect  the  com¬ 
mon  species  flow  field  may  have  a  pronounced  affect 
on  the  raire  species  flow  field. 

In  order  to  achieve  this  coupling,  the  particle 
properties  are  binned  in  every  cell  during  the  sam¬ 
pling  period  of  the  base  simulation.  Five  different 
properties  are  binned:  the  three  components  of  ve¬ 
locity,  the  rotational  energy  and  the  vibrational  en¬ 
ergy.  During  the  overlay  simulation,  a  temporary 
common  particle  is  created  to  act  as  a  collision  part¬ 
ner  for  each  rare  particle  in  every  cell.  These  tem¬ 
porary  particles  have  their  properties  determined 
through  sampling  of  the  bins  stored  during  the  base 
simulation.  From  these  pairs,  the  appropriate  num¬ 
ber  of  common- rare  collisions  is  determined.  After 
standard  DSMC  collision  mechanics  is  performed, 
the  temporary  particles  are  discarded.  Thus,  the 
particles  of  the  overlay  simulation  interact  with  the 
background  flow  field  at  the  microscopic  level.  It 
is  important  to  note  that  given  the  assumption  that 
the  rare  species  have  a  number  density  that  is  several 
orders  of  magnitude  lower  than  the  common  species, 
the  probability  of  common-rare  collisions  is  much 
greater  than  rare-rare  collisions.  In  the  cases  con¬ 
sidered  here,  rare-rare  collisions  are  neglected  in  the 
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overlay  simulation. 

Chemistry  in  DSMC  Overlay 

Applying  the  DSMC  overlay  procedure  outlined 
above  to  chemically  reacting  flows  is  somewhat  more 
involved.  It  is  important  to  determine  what  classes 
of  collisions  will  lead  to  chemical  reactions  that  affect 
the  rare  species.  The  three  possible  candidates  are 
common-common  collisions,  common-rare  collisions 
and  rare-rare  collisions.  As  was  previously  men¬ 
tioned,  rare-rare  collision  are  neglected,  so  no  reac¬ 
tions  stemming  from  these  types  of  collision  are  con¬ 
sidered.  Common- rare  collisions  are  handled  using 
the  bins  and  temporary  particles.  Standard  DSMC 
reaction  routines  are  used  to  handle  the  probabil¬ 
ity  and  mechanics  of  common-rare  reactions.  How¬ 
ever,  consideration  must  be  given  to  the  relative  im¬ 
portance  of  common-rare  reactions  compared  to  the 
common-common  reactions.  Since  the  probability 
of  a  common-rare  collision  is  several  orders  of  mag¬ 
nitude  lower  than  a  common-common  collision,  the 
reaction  probability  must  be  very  high  in  order  for 
the  common-rare  reaction  to  be  important  relative 
to  the  common-common  reaction.  In  the  cases  con¬ 
sidered  here,  the  reaction  probabilities  for  the  reac¬ 
tions  of  interest  are  not  sufficiently  high,  and  the 
inclusion  of  common-rare  reaction  has  no  quantita¬ 
tive  effect  on  the  rare  species  flow  field.  This  leaves 
only  common-common  reactions  as  the  important 
reactions  in  terms  of  the  production  of  rare  species. 

Common-common  reactions  cannot  be  handled 
during  the  overlay  simulation  using  standard  DSMC 
techniques  as  no  common-common  collision  pairs  are 
processed.  Regardless  of  the  technique  used  to  de¬ 
termine  the  concentration  and  properties  of  the  rare 
particles  created  via  these  reactions,  their  introduc¬ 
tion  into  the  overlay  simulation  will  take  the  form  of 
a  source  term  for  the  different  rare  species  in  every 
cell.  A  completely  general  scheme  for  accomplishing 
this  would  require  binning  more  detailed  information 
during  the  base  simulation.  For'  the  work  presented 
here,  a  more  approximate  method  of  determining 
both  the  rate  of  production  and  the  properties  of 
the  rare  species  due  to  common-common  reactions 
is  used.  The  predominant  reaction  considered  is: 

N2  +  O  NO  -h  N 

It  is  assumed  that  only  the  most  energetic  collisions 
lead  to  chemical  reaction.  In  the  Atmospheric  Ex¬ 
plorer  flows,  this  occurs  when  a  free  stream  particle 
collides  with  a  particle  that  has  rebounded  from  the 
ram  surface  of  the  satellite.  Thus,  there  are  two  dif¬ 
ferent  types  of  collisions  leading  to  the  formation  of 


nitric  oxide: 

N2  +  0  --4  nT)  +  N  (I) 

N2  +  O  NO  -h  N  (II) 

where  the  arrows  indicate  the  direction  in  which 
the  molecule  is  moving.  Thus,  reaction  I  represents 
a  free  stream  nitrogen  molecule,  moving  at  orbital 
speed,  colliding  with  an  atomic  oxygen  molecule  that 
has  reflected  from  the  surface.  Reaction  II  represents 
free  stream  atomic  oxygen  reacting  with  reflected 
N2. 

Data  on  the  velocity  distributions  of  surface  scat¬ 
tered  8  km/s  0  and  N2  is  limited.  Examination  of 
the  collision  of  8  km/s  oxygen  atoms  with  engineer¬ 
ing  surfaces  indicates  that  the  scattering  is  largely 
diffuse,  although  the  scattered  atoms  have  a  broad 
velocity  distribution  centered  around  1  km/s}^  Lim¬ 
ited  measurements  of  high  velocity  N2  scattering 
suggest  that  the  scattered  velocity  is  approximately 
2  km/s.^^  Using  a  value  of  8  km/s  for  orbital  speed, 
this  results  in  a  relative  velocity  of  9  km/s  or  10  km/s 
for  Reactions  I  and  II,  respectively.  Using  these  ap¬ 
proximations  and  appropriate  values  for  the  cross- 
section  of  the  reaction  at  these  speeds,  the  rate  of 
production  of  NO  and  N  can  be  determined  in  each 
cell  using  the  equation: 

dn(NO)  dn(N) 

Note  that  the  value  for  the  number  density  used  for 
nitrogen  and  atomic  oxygen  must  be  modified  to  ac¬ 
count  for  the  assumptions  that  lead  to  the  choice  of 
values  for  the  relative  velocity.  An  appropriate  vadue 
for  a  is  found  from  consideration  of  other  sources. 

In  each  cell  at  each  time  step,  a  number  of  NO 
and  N  particles  are  generated  according  to  the  pro¬ 
duction  term  evaluated  using  number  densities  from 
the  bcLse  simulation.  For  the  non-flow  direction  com¬ 
ponents  of  velocity,  the  properties  of  the  NO  and  N 
particles  created  can  be  approximated  from  the  ve¬ 
locity  distributions  of  N2  and  O,  respectively.  The 
velocity  component  in  the  direction  of  the  flow  for 
the  newly  formed  NO  and  N  is  approximated  by  the 
center  of  mass  velocity  of  the  collision,  with  a  ther¬ 
mal  spread  obtained  by  sampling  the  azimuthal  dis¬ 
tribution  in  that  cell. 

Results 

Simulations  are  performed  at  10  km  intervals  over 
the  altitude  range  of  140  km  to  180  km  for  flow 
over  the  Atmospheric  Explorer  satellite.  The  ge¬ 
ometry  of  the  satellite  is  a  0.7  m  radius  cylinder, 
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h  [km] 

Kn 

n(NO)  [m-^] 

n(N)  [m 

140 

27.5 

7.0el2 

7.54el2 

150 

47.7 

5.5el2 

1.39el3 

160 

75.9 

4.0el2 

2.17el3 

170 

114.0 

3.0el2 

2.90el3 

180 

164.5 

2.0el2 

3.44el4 

Table  1:  Summary  of  Free  Stream  Condition 


1.0  ra  in  length.  At  each  altitude,  three  different 
simulations  are  made.  The  first  of  these  is  a  base 
DSMC  simulation,  where  the  flow  field  of  the  com¬ 
mon  species  is  determined  and  particle  properties 
are  binned.  Two  overlay  DSMC  simulations  axe  then 
run  using  the  base  simulation  for  the  background 
conditions.  These  simulations  differ  in  the  choice 
of  cross-sections  for  the  source  term  chemical  reac¬ 
tions  for  NO  and  N.  In  the  first  overlay  simulation, 
cross-sections  of  cr/  =  4  x  10“*^^  m^  for  Reaction  I 
and  an  =  6  x  m^  for  Reaction  II  are  used.^® 

In  the  other  overlay  simulation,  values  for  aj  and 
an  obtained  from  a  molecular  dynamics  trajectory 
analysis  are  used^^  and  are  five  times  greater  than 
the  values  used  in  the  first  overlay  simulation.  A 
summary  of  free  stream  values  is  given  in  Table  1. 


Computational  Performance 

The  simulations  presented  here  are  performed  us¬ 
ing  a  highly  modified  version  of  MONACO,  a  paral¬ 
lel,  object-oriented  DSMC  code  developed  at  Cornell 
University.^®  The  code  is  run  on  eight  nodes  of  an 
IBM  SP-2.  The  total  number  of  particles  used  dur¬ 
ing  the  simulation  varies  from  150,000  to  1,000,000, 
depending  on  the  altitude.  For  the  base  simula¬ 
tions,  the  flow  takes  6,000  time  steps  to  reach  steady 
state  and  the  flow  field  is  sampled  for  2,000  time 
steps.  The  transient  period  for  overlay  simulations 
also  takes  6,000  times  steps.  For  overlay  simulations, 
5,000  sampling  steps  are  used  in  order  to  accurately 
resolve  wall  distributions.  The  computational  effort 
required  also  varies  significantly  between  the  vari¬ 
ous  runs.  Base  simulations  are  runs  with  a  compu¬ 
tational  cost  of  10  /is/particle/time  step.  The  ad¬ 
ditional  overhead  of  sampling  the  velocity  distribu¬ 
tions  increases  the  cost  to  24  ^s/particle/time  step 
in  the  overlay  simulation.  This  variation  in  com¬ 
putational  cost  coupled  with  the  varying  number  of 
particles  results  in  a  range  of  CPU  times  from  2  to 
10  hours. 


Figure  1:  Computational  Grid  Used  At  All  Altitudes 

General  Flow  Field  Structure 

A  grid  consisting  of  722  computational  cells  is  used 
for  all  altitudes  and  is  shown  in  Figure  1.  Cells  of  a 
typical  DSMC  grid  scale  in  size  with  the  mean  free 
path  of  the  flow.  For  the  altitudes  considered  the 
mean  free  path  varies  from  about  3  m  to  200  m  and 
thus  is  not  a  useful  scaling  parameter.  Instead,  cell 
sizes  are  scaled  in  order  to  achieve  useful  resolutions 
during  sampling.  This  permits  the  same  grid  to  be 
used  for  all  the  altitudes  and  cases  considered. 

The  overall  structure  of  the  flow  field  is  given  in 
Figure  2,  which  shows  the  translational  temperature 
at  140  km.  A  diffuse  bow  shock  can  be  seen  in  front 
of  the  body  of  the  Atmospheric  Explorer,  consistent 
with  rarefied,  hypersonic  flow.  No  quantitative  con¬ 
clusions  should  be  drawn  from  plots  such  as  this, 
as  the  relevance  of  thermodynamic  variables  in  this 
rarefied  regime  is  questionable.  It  is  provided  more 
for  a  qualitative  overview  of  the  structure  of  the  flow 
field. 

Plots  of  the  number  density  for  all  five  species 
along  the  axis  at  140  km  and  using  the  first  set  of 
cross-sections  are  shown  in  Figure  3.  The  data  for 
these  density  profiles  is  obtained  from  both  the  base 
and  overlay  simulations.  Immediately  clear  is  the 
large  difference  in  the  number  density  between  the 
common  species  and  the  rare  species. 

Sensitivity  to  Altitude 

Effects  of  varying  the  altitude  on  the  form  of  the  ve¬ 
locity  distributions  is  now  considered.  For  the  distri¬ 
butions  presented,  the  comparison  is  made  between 
cases  at  140  km  and  180  km,  using  the  first  cross- 
section.  Above  this  altitude  range,  the  concentra¬ 
tion  of  atomic  nitrogen  has  increased  to  such  a  level 


5 

American  Institute  of  Aeronautics  and  Astronautics 


AlAA  96-1847 


Axial  Distance  (m) 

Figure  2:  Common  Species  TransIationcJ  Tempera¬ 
ture  at  140  km 


I  .I . . . 1  .  .  .  >  i 

'^-10  -9  -8  -7  -6  -5  -4  -3  -2  -1 


Axial  Distance  (m) 

Figure  3:  Number  Density  Profiles  Along  the  Axis 
at  140  km 

that  it  can  no  longer  be  considered  a  rare  species. 
Due  to  the  affects  this  would  have  on  the  source 
term  chemistry,  simulations  above  this  altitude  were 
not  performed.  Unless  otherwise  stated,  the  distri¬ 
butions  from  the  intermediate  altitudes  follow  the 
trends  from  lower  to  higher  altitude  described  in  the 
plots.  Note  that  in  all  cases,  the  distributions  are 
normalized.  Thus,  an  increase  in  one  area  of  the 
distributions  must  lead  to  a  corresponding  decrease 
in  other  areas. 

Figure  4  illustrates  the  normal  velocity  distribu¬ 
tions  of  nitric  oxide  particles  impacting  the  surface. 
At  180  km,  the  majority  of  the  particles  striking  the 
surface  have  a  velocity  centered  around  8  kra/s  with 
a  thermal  spread  corresponding  to  free  stream  con¬ 
ditions.  This  peak  represents  free  stream  patrticles 
that  have  reached  the  surface  of  the  satellite  with- 
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Figure  4:  Variation  of  the  Normal  Velocity  Distri¬ 
bution  Functions  at  the  Wall  with  Altitude 

out  undergoing  any  collisions.  The  two  other  peaks 
in  the  180  km  distribution  are  representative  of  the 
nitric  oxide  particles  created  from  source  term  reac¬ 
tions.  The  two  distinct  peaks  are  due  to  the  different 
center  of  mass  velocities  of  the  two  production  reac¬ 
tions.  The  velocity  distribution  at  140  km  shows  a 
similar  form.  The  main  distinction  lies  in  the  rela¬ 
tive  importance  of  the  source  term  NO  to  free  stream 
NO.  At  140  km,  the  majority  of  NO  particles  striking 
the  surface  are  created  in  the  flow  field,  as  opposed 
to  those  present  in  the  ambient  atmosphere. 

Some  further  comments  about  the  normal  velocity 
distributions  are  warranted.  The  distributions  have 
not  been  plotted  for  the  lower  range  of  velocities 
(vn  <  3  km/s)  because  the  distributions  are  iden¬ 
tically  zero  in  this  range.  This  indicates  that  the 
number  of  rare-common  collisions  in  the  flow  field  is 
extremely  limited,  even  at  the  lowest  of  the  altitudes 
considered.  If  there  were  a  substantial  number  of 
rare-common  collisions,  a  discernible  number  of  NO 
particles  striking  the  surface  would  have  a  low  nor¬ 
mal  velocity.  It  should  also  be  noted  that,  due  to 
the  importance  of  the  source  term  chemistry  at  the 
lower  altitudes,  the  average  velocity  of  the  NO  par¬ 
ticles  hitting  the  surface  is  comparatively  lower  than 
at  the  higher  altitudes.  This  is  important  in  inter¬ 
preting  the  effects  of  using  the  larger  cross-sections. 

Sensitivity  to  Reaction  Cross-Sections 

Comparison  of  the  effects  of  the  two  reaction  cross- 
sections  tested  is  most  appropriately  studied  at  the 
lowest  altitude  considered,  140  km.  At  this  altitude 
the  base  flow  field  density  is  the  highest,  and  there¬ 
fore  the  effects  of  changing  the  cross-section  are  the 
most  noticeable.  Figures  5-8  compare  the  nitric 
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Figure  5:  Variation  of  the  X- Velocity  Distribution 
Functions  with  Cross-section  in  an  Inflow  Cell 

oxide  X- velocity  distributions  for  cases  using  the  two 
cross-sections  at  140  km  at  four  points  in  the  flow 
field.  The  four  points  chosen  are  the  inflow  cell  along 
the  axis  approximately  10  m  from  the  the  front  of 
the  satellite,  a  mid-flow  field  cell  approximately  3  m 
from  the  surface  on  the  axis,  the  last  cell  before  the 
wail  on  the  axis  and  on  the  ram  surface  of  the  satel¬ 
lite  itself. 

Figures  5  compares  the  x-velocity  distribution  in 
the  inflow  cell  along  the  axis.  In  both  cases,  two 
distinct  groups  of  particles  axe  evident.  The  group 
with  the  high,  positive  x-velocity  represents  free 
stream  particles  moving  towards  the  satellite  at  or¬ 
bital  speeds.  The  other  represents  particles  that 
have  collided  with  the  front  face  of  the  satellite 
and  are  moving  away  with  a  negative  x-velocity. 
Given  the  extremely  rarefied  conditions,  these  par¬ 
ticles  travel  throughout  the  entire  flow  field  with¬ 
out  undergoing  collisions.  Note  that  the  relative  fre¬ 
quency  of  particles  in  the  rebounded  group  is  larger 
with  the  second,  larger  cross-sections  ((72)-  With 
the  larger  cross-sections,  more  particles  have  been 
created  throughout  the  flow  field  due  to  the  source 
chemistry.  These  particles  also  collide  with  the  sur¬ 
face,  rebound  and  are  sampled  in  the  inflow  cell  with 
a  negative  x-velocity. 

A  comparison  of  the  x-velocity  distributions  in  the 
mid-flow  field  cell  on  the  axis  is  presented  in  Fig¬ 
ure  6.  The  two  groups  evident  in  the  inflow  cell  are 
also  seen  here.  The  relative  increase  of  the  negative 
x-velocity  group  with  the  second  cross-section  is  also 
apparent.  In  addition  to  these  groups,  two  smaller 
groups  are  distinguishable.  Noticeably  smaller  in 
magnitude  than  the  first  two  groups,  these  can  be 
seen  at  an  x-velocity  of  about  5  km/s.  These  groups 
represent  the  particles  created  through  the  source 


Figure  6:  Variation  of  the  X- Velocity  Distribution 
Functions  with  Cross-section  in  an  Mid-Flow  Field 
Ceil 

term  chemistry  that  have  yet  to  collide  with  the  AE 
satellite.  The  two  peaks  in  this  grouping  show  the 
difference  in  the  center  of  mass  velocity  in  the  two 
source  term  reactions.  When  an  N2  molecule  travel¬ 
ing  at  orbital  speeds  collides  with  a  reflected  oxygen 
atom,  the  center  of  mass  velocity  is  larger  than  when 
the  atomic  oxygen  is  traveling  at  orbital  speeds. 

Near  the  wall,  the  comparison  between  the  dis¬ 
tributions  is  noticeably  different.  In  both  cases 
shown  in  Figure  7,  the  distribution  is  dominated  by 
a  group  of  reflected  particles,  moving  in  the  negative 
x-direction  with  a  relatively  slow  speed.  The  form 
of  these  distributions  is  consistent  with  the  number 
density  profiles.  The  magnitude  of  the  velocity  of 
the  rebounding  particles  is  relatively  small,  so  that 
it  will  take  a  much  longer  time  for  these  particles 
to  leave  the  cell  than  it  will  for  the  high  speed  free 
stream  particles  in  the  opposite  direction.  Conse¬ 
quently,  at  any  time  there  are  many  more  reflected 
particles  than  free  stream  particles  in  cells  near  the 
surface. 

The  final  comparison  of  the  effects  of  changing 
the  cross-section  is  given  in  Figure  8,  which  shows 
normal  velocity  distributions  along  the  wall.  In  both 
distributions,  the  presence  of  source  term  NO  as  well 
as  free  stream  NO  is  evident.  Also  distinguishable 
is  the  bimodal  behavior  of  the  source  term  parti¬ 
cles.  Once  again,  the  increased  importance  of  the 
source  term  chemistry  is  visible;  the  distribution  in 
the  second  cross-section  case  is  heavily  weighted  to¬ 
wards  these  particles.  Note  that  the  predominance 
of  source  term  particles  is  so  large  that  it  greatly 
overshadows  the  free  stream  peak. 
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Figure  7:  Variation  of  the  X- Velocity  Distribution 
Functions  with  Cross-section  in  a  Near  Wall  Cell 


Figure  8:  Variation  of  the  Normal  Velocity  Distri¬ 
bution  Functions  with  Cross-section  at  the  Wall 

Glow  Brightness  Predictions 

Figures  9  and  10  show  a  comparison  between  exper¬ 
imental  AE  data  for  the  glow  brightness  at  6563  A 
at  varying  altitudes^  and  predictions  based  on  the 
DSMC  overlay  simulations  using  the  cross-sections 
sets  cji  and  02^  In  both  cases,  the  results  have  been 
scaled  to  units  of  Rayleighs,  where  1  Rayleigh  is 
defined  as  10®  photons/cm^/s.  The  DSMC  over¬ 
lay  prediction  of  the  glow  brightness  is  based  on  the 
average  nitric  oxide  flux  hitting  the  surface.  It  is  as¬ 
sumed  that  each  NO  molecule  striking  the  surface  in¬ 
teracts  with  an  adsorbed  oxygen  atom,  forms  a  NO2 
molecule  in  the  gas  phase  that  spontaneously  decays, 
emitting  a  photon.  The  NO2  molecule  is  assumed 
to  emit  the  photon  in  a  random  direction,  with  all 
directions  being  equally  probable.  The  DSMC  pre¬ 


diction  of  the  brightness  is  further  scaled  so  that  the 
prediction  agrees  with  experiment  at  an  altitude  of 
160  km.  The  magnitude  of  this  scaling  is  substan¬ 
tial,  on  the  order  of  lO'*,  and  is  larger  for  the  second 
cross-section  cases.  This  scaling  is  required  for  the 
most  part  to  correct  for  the  unit  probability  assump¬ 
tion  for  the  gas-surface  interaction.  A  true  reaction 
probability  is  likely  to  be  significantly  less  than  one, 
so  that  the  true  production  rate  of  NOJ  is  much  less 
than  the  incident  NO  flux  striking  the  surface.  In 
each  case,  a  prediction  of  the  glow  brightness  based 
solely  on  free  stream  nitric  oxide  is  included  for  com¬ 
parison.  At  a  particular  altitude,  the  magnitude  of 
the  difference  between  the  DSMC  overlay  prediction 
and  the  free  stream  estimate  is  an  indication  of  the 
importance  of  gas-gas  collisions  and  chemistry. 

DSMC  predictions  for  the  first  cross-sections  are 
presented  in  Figure  9.  In  general,  the  predictions 
show  good  agreement  over  the  range  of  altitudes 
considered.  The  deviation  from  the  free  stream  ref¬ 
erence  is  greatest  at  the  lowest  altitude.  This  is 
consistent  with  the  input  conditions  and  the  pre¬ 
vious  results.  At  the  lower  altitudes,  the  density  is 
higher  so  the  production  from  source  term  chemistry 
is  greater.  This  results  in  a  higher  overall  NO  flux 
which  in  turn  leads  to  increased  glow  production.  As 
altitude  increases,  the  DSMC  predictions  approach 
the  free  stream  reference  predictions.  At  the  high¬ 
est  altitude  considered,  180  km,  the  two  predictions 
almost  coincide.  This  is  an  indication  that  produc¬ 
tion  of  nitric  oxide  in  the  flow  field  is  unimportant 
at  greater  altitudes  and  that  the  glow  is  almost  en¬ 
tirely  due  to  free  stream  nitric  oxide.  In  addition,  it 
demonstrates  the  importance  of  the  initial  estimates 
of  the  nitric  oxide  densities,  as  the  DSMC  predic¬ 
tions  will  inevitably  match  the  free  stream  predic¬ 
tions  at  higher  altitudes.  Error  in  the  nitric  oxide 
number  density  directly  translates  into  error  in  the 
glow  production. 

Figure  10  gives  the  predictions  based  on  the  sec¬ 
ond  cross-sections.  The  DSMC  simulations  are  seen 
to  overpredict  the  brightness  at  the  lower  altitudes, 
where  the  source  term  production  of  NO  is  more  im¬ 
portant.  Larger -predicted  brightness  using  the  sec¬ 
ond  cross-sections  is  certainly  to  be  expected,  as  the 
larger  magnitude  of  the  second  cross-sections  results 
in  an  increase  in  the  production  of  NO.  It  is  not 
entirely  clear,  however,  that  the  first  cross-section 
should  be  preferred  over  the  second  simply  on  the 
basis  of  the  superior  agreement  with  the  experimen¬ 
tal  data.  The  effects  of  the  normal  velocity  distri¬ 
butions  along  the  wall  for  the  incident  nitric  oxide 
is  neglected  in  the  simulations  presented  here.  For 
the  lower  altitude  case,  as  seen  in  Figure  4,  a  larger 
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Figure  9;  Comparison  of  Predicted  Glow  Brightness 
Using  (Ti  with  AE  Experimental  Data 


Figure  10:  Comparison  of  Predicted  Glow  Bright¬ 
ness  Using  ct2  with  AE  Experimental  Data 

fraction  of  the  NO  striking  the  surface  is  produced 
through  the  source  term  chemistry.  This  results  in 
the  NO  striking  the  surface  having  a  lower  average 
velocity  at  the  lower  altitudes.  If  the  simple  assump¬ 
tion  is  made  that  the  higher  the  incident  energy,  the 
higher  the  probability  of  reaction,  then  the  predicted 
brightness  at  the  lower  altitudes  would  be  reduced. 
So,  taking  into  account  the  effect  of  the  normal  ve¬ 
locity  distribution  may  result  in  the  second  cross- 
section  giving  better  agreement  with  experiment. 

Conclusions 

A  DSMC  overlay  technique  suitable  for  simulating 
rare  species  has  been  developed.  This  method  was 
applied  to  the  flow  around  the  Atmospheric  Explorer 


satellite  in  order  to  predict  the  NO2  glow  brightness. 

Feasibility  of  the  DSMC  overlay  technique  has 
been  demonstrated.  Recall  the  number  density  pro¬ 
files  shown  in  Figure  3  indicate  that  the  relative  dif¬ 
ference  between  the  densities  of  the  common  and 
rare  species  is  three  orders  of  magnitude.  Thus,  the 
DSMC  overlay  method  can  clearly  resolve  the  behav¬ 
ior  of  rare  species.  Numerous  velocity  distributions 
for  nitric  oxide  have  also  been  presented,  showing  the 
ability  of  the  method  to  maintain  microscopic  detail 
for  the  rare  species  throughout  the  simulation. 

Prediction  of  the  glow  brightness  c^er  a  range 
of  altitudes  has  been  made.  For  the  reaction 
cross-sections  based  on  measurement,  the  qualitative 
agreement  of  the  scaled  predictions  with  the  exper¬ 
iment  data  was  quite  good.  For  the  cross-sections 
based  on  computational  chemistry,  the  agreement 
was  not  as  good  at  the  lower  altitudes  where  the 
prediction  was  significantly  larger  than  experiment. 
This  overprediction  could  be  tempered  by  the  inclu¬ 
sion  of  collisional  energy  effects  in  the  formation  of 
NO2.  In  both  cases,  the  importance  of  source  term 
chemist rv  was  seen  to  diminish  at  higher  altitudes. 
At  adtitudes  greater  than  200  km,  the  results  suggest 
that  the  glow  production  is  entire  due  to  ambient 
nitric  oxide.  This  highlights  the  importance  of  ac¬ 
curately  determining  the  ambient  concentrations  of 
the  rare  species. 
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Abstract 

The  Skipper  flight  experiment  is  due  for  launch  in 
December  1995.  The  experiment  is  designed  to  make 
measurements  of  visible  and  ultra-violet  emissions 
caused  by  the  interaction  of  a  hypersonic  vehicle 
with  the  atmosphere  at  very  high  altitudes.  The 
mission  consists  of  two  distinct  phases:  (1)  measure¬ 
ment  of  spacecraft  glow  emission  during  orbit  of  the 
earth  over  the  altitude  range  of  300  down  to  about 
130  km:  and  (2)  measurement  of  bow-shock  gener¬ 
ated  radiation  during  reentry  over  the  altitude  range 
of  180  to  about  100  km.  This  document  provides 
pre- flight  predictions  for  emission  generated  in  both 
phases  of  the  flight.  The  purpose  of  these  predic¬ 
tions  is  to  assist  during  the  execution  of  the  Skipper 
flight-  The  results  should  be  regarded  as  prelimi¬ 
nary.  It  is  anticipated  that  more  detailed  modeling 
will  be  performed  after  the  flight  based  on  the  data 
provided  by  Skipper. 

Introduction 

The  Skipper  flight  experiment  is  designed  to  char¬ 
acterize  the  aerothermochemistry  experienced  by  a 
spacecraft  traveling  at  hypersonic  speeds  during  low 
earth  orbit  and  atmospheric  reentry.  The  flight  is 
composed  of  two  distinct  stages.  The  first  stage 
consists  of  a  number  of  elliptical,  low-earth  orbits, 
during  which  spacecraft  glow  measurements  will  be 
taken.  In  the  second  stage,  the  satellite  will  take 
radiation  measurements  during  reentry. 


Earlier  flight  experiments  have  successfully  gath¬ 
ered  data  for  use  in  refining  air-chemistry  modeling. 
One  of  these  experiments  (BSU\'-1)  operated  at  a 
velocity  of  about  3.5  km/s  and  covered  an  altitude 
range  of  40  km  to  70  km.“  While  initial  calculations 
of  the  ultraviolet  emissions  were  low  by  several  or¬ 
ders  of  magnitude,  improvements  to  the  air  radiation 
models  give  good  agreement  between  the  measured 
and  calculated  emissions,  A  second  flight  experi¬ 
ment  (BSUV-2)  operated  at  a  velocity  of  5.1  km/s 
over  the  altitude  range  of  65  km  to  100  km.'^  Due  to 
the  rarefied  nature  of  the  atmosphere  at  the  higher 
altitudes,  direct  simulation  Monte  Carlo  (DSMC) 
calculations  give  better  agreement  with  measured 
emissions  of  XO  than  CFD  predictions.  Under  these 
conditions,  simulation  is  required  both  of  species 
present  in  extremely  low  concentrations  and  of  ex¬ 
tremely  rare  collision  events.  Initial  DSMC  calcula¬ 
tions  lacking  the  particle  resolution  to  simulate  such 
rare  species  and  events  grossly  underpredicted  the 
X^O  emission.  Use  of  an  overlay  technique  in  con¬ 
junction  with  a  standard  DSMC  calculation  allowed 
both  common  and  rare  species  and  events  to  be  mod¬ 
eled  correctly,  bringing  the  predictions  much  closer 
to  the  measured  values.’^ 

This  report  presents  preflight  predictions  of  ra¬ 
diation  measurements.  These  predictions  are  not 
meant  to  be  definitive,  but  rather  are  intended  for 
assistance  during  the  operation  of  the  Skipper  space¬ 
craft.  It  is  anticipated  that  the  predictions  discussed 
here  will  give  only  qualitative  agreement  with  the  ac¬ 
tual  data  set.  Predictions  are  provided  for  both  the 
orbital  and  reentry  portions  of  the  Skipper  mission. 
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The  results  and  the  methods  employed  are  presented 
in  the  following  sections. 

UV  and  VUV  Emissions  Ob¬ 
served  During  Reentry 

In  this  section  predictions  are  provided  of  ultra¬ 
violet  and  vacuum-ultra-violet  radiative  emissions 
for  nitric  oxide  (NO),  and  atomic  oxygen  (0).  Re¬ 
sults  are  presented  over  the  altitude  range  from  160 
down  to  100  km.  We  describe  the  solution  procedure 
for  computation  of  the  flow  fields  around  the  Skip¬ 
per  vehicle,  the  chemistry  models,  and  the  radiation 
models. 

Solution  Procedure 

A  multistep  procedure  is  used  to  calculate  rtidiative 
emission  predictions,  consisting  of  a  DSMC  calcula¬ 
tion,  an  overlay  calculation,  and  a  radiance  calcu¬ 
lation.  The  DSMC  simulation  takes  as  input  free- 
stream  conditions  corresponding  to  the  altitude  in 
question.  Due  to  limitations  in  particle  resolution, 
this  calculation  only  simulates  the  major  species. 
N2,  O2  and  O,  to  determine  macroscopic  proper¬ 
ties  such  as  temperature  and  mass  density  in  the 
flow  field.  The  resulting  solution  serves  as  input  to 
the  overlay  code.  The  overlay  calculation  determines 
the  concentrations  of  the  trace  species,  NO  and  N. 
As  with  previous  use  of  the  overlay  technique,  the 
trace  species  are  assumed  to  be  present  in  such  low 
concentrations  that  reactions  involving  these  species 
have  no  effect  on  the  overall  temperature  and  den¬ 
sity  of  the  flow  field,  nor  on  the  concentrations  of  the 
major  species.  Finally,  the  complete  flow  field  solu¬ 
tion,  including  both  major  and  trace  species,  is  used 
in  the  NEQAIR  code^  ®  to  compute  the  radiative 
emissions. 

Atmospheric  Model 

Two  separate  atmospheric  models  are  used.  The 
first  uses  the  MSISE-90  model'  to  predict  the  at¬ 
mospheric  structure.  Inputs  to  the  model  include 
date,  time,  latitude,  longitude,  solar  flux  indices  and 
the  geomagnetic  index.  Predictions  of  atmospheric 
structure  are  obtained  for  latitude  23°N  and  lon¬ 
gitude  164°E,  an  approximate  location  for  reentry. 
The  second  atmosphere  used  is  simply  a  standard  at¬ 
mosphere  taken  from  the  Handbook  of  Geophysics 
and  the  Space  Environment.^  While  not  specific  to 
the  reentry  point  or  any  of  the  other  parameters  that 
determine  atmospheric  composition,  the  standard 


atmosphere  does  form  a  good  baseline  to  compare 
with  the  MSIS  atmosphere. 

Determination  of  the  concentrations  of  N  and  NO 
is  somewhat  more  complicated.  While  the  MSIS 
model  gives  variation  of  N  number  density  with  al¬ 
titude.  the  standard  atmosphere  does  not  include 
trace  species.  The  MSIS  value  for  the  number  den¬ 
sity  of  N  is  used  as  the  standard  atmosphere  value  as 
well.  Because  the  concentration  of  the  trace  species 
is  much  smaller  than  the  major  species,  estimating 
the  N  concentration  in  this  manner  does  not  affect 
the  DSMC  simulation.  The  concentration  of  NO  is 
so  low  compared  to  the  major  species  that  the  MSIS 
atmosphere  does  not  predict  NO  number  density. 
Nitric  oxide  number  density  for  both  atmospheres 
is  estimated  from  data  presented  by  McCoy.^  At  all 
altitudes  in  question,  the  amount  of  NO  produced 
through  chemical  reaction  far  outweighs  the  amount 
present  in  the  free-stream.  Therefore,  reasonable  er¬ 
rors  in  these  estimates  should  not  affect  the  emission 
predictions. 

DSMC  Calculations 

The  code  employed  for  the  DSMC  calculations  is 
based  on  that  described  by  Boyd  and  Gokcen.^°  The 
satellite  is  assumed  to  be  roughly  cylindrical,  with  a 
diameter  of  O.S  m  and  an  axial  length  of  1.25  in.  For 
all  simulations,  the  base  of  the  cylinder  is  positioned 
at  the  origin,  with  the  body  .set  in  the  negative  x- 
direction.  Due  to  the  extreme  change  in  atmospheric 
density  over  the  altitude  range  in  question,  a  differ¬ 
ent  grid  is  used  at  each  altitude.  Grid  sizes  range 
from  89  x  36  cells,  covering  an  actual  domain  of  ap¬ 
proximately  2  in  X  1  in.  for  the  100  km  case,  to 
19  X  45  cells,  covering  a  physical  domain  of  approx¬ 
imate!  v  20  in  x  15  in.  for  the  160  km  case.  Due  to 
the  large  difference  in  the  concentrations  of  the  ma¬ 
jor  species  (NT,  02,0)  to  the  trace  species  (NO.N), 
which  varies  from  4  to  6  orders  of  magnitude  in  the 
cases  considered,  only  reactions  involving  the  ma¬ 
jor  species  are  simulated  in  the  DSMC  code.  Out¬ 
put  of  the  DSMC  simulation  gives  concentrations 
for  the  major  species,  as  well  as  temperatures  and 
mass  density  values  for  each  cell.  While  reactions 
involving  the  trace  species  are  not  simulated  due  to 
lack  of  particle  resolution,  it  is  assumed  that  these 
species  are  indeed  rare  enough  that  including  them 
later  in  the  overlay  calculation  will  not  change  the 
macroscopic  flow  field  values. 

Overlay  Calculations 

Following  the  method  that  was  successful  in  simu¬ 
lating  trace  species  reactions  for  the  BSUV-2  flight,*^ 
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Figure  1:  N2  and  Atomic  O  Concentration  with  Al¬ 
titude  for  the  MSIS  and  Standard  Atmospheres 


an  overlay  code  is  employed  using  the  output  flow 
field  from  the  DSMC  simulation.  The  overlay  code 
was  written  by  Candler,  with  minor  modifications 
made  for  the  present  study.  As  mentioned  above, 
the  overlay  method  assumes  that  the  trace  species 
are  present  in  such  small  amounts  that  they  do  not 
affect  the  overall  flow  field.  The  DSMC  calculated 
flow  field  serves  as  the  unvarying  flow  field.  Mass 
conservation  is  then  solved  using  the  following  con¬ 
servation  equation: 


for  each  of  the  n  trace  species,  is  the  trace 
species  density,  Vsj  is  the  diffusion  velocity  and 
is  the  chemical  source  term  for  species  .s.  The  mass- 
averaged  velocity,  Uj,  is  determined  from  the  DSMC 
solution. 


Chemical  Reactions 

One  of  the  primary  measurements  on  the  Skipper 
flight  is  the  230  nm  ultraviolet  emission  of  nitric  ox¬ 
ide.  Accordingly,  it  is  possible  that  the  modeling 
of  any  reactions  resulting  in  the  production  of  NO 
could  have  a  great  effect  on  the  NO  emission  pre¬ 
dictions,  For  rarefied  air.  the  two  most  important 
NO  formation  reactions  are  the  Zeldovich  exchange 
reactions: 


of  No  is  extremely  rare  and  N  is  present  only  as  a 
trace  species  in  the  free-stream.  so  there  is  gener¬ 
ally  very  little  atomic  N  available  to  drive  this  re¬ 
action.  Therefore,  the  second  of  the  listed  exchange 
reactions  is  much  more  important  in  controlling  the 
production  of  NO.  While  the  first  exchange  reac¬ 
tion  is  modeled,  the  corresponding  rate  coefficient 
is  not  varied  in  this  study.  In  the  second  exchange 
reaction,  the  concentration  of  atomic  0  controls  the 
production  of  NO.  In  the  BSUV-2  case,  it  was  found 
that  the  dissociation  rate  of  O2  had  to  be  mod¬ 
eled  very  carefully  in  order  to  accurately  predict 
the  amount  of  NO  produced.  In  the  higher  alti¬ 
tude  range  being  considered  here,  dissociation  mod¬ 
eling  is  less  important  as  there  is  already  a  significant 
amount  of  atomic  O  present  in  the  ambient  atmo¬ 
sphere:  solar  radiation  results  in  a  large  amount  of 
photodissociation  of  Oo.  It  is  therefore  important 
that  the  forward  rate  of  the  second  exchange  reac¬ 
tion  be  modeled  accurately.  Two  different  reaction 
rates  are  used  in  making  the  predictions.  The  first  of 
these  is  the  rate  proposed  by  Park.^^  which  assumes 
an  Arrhenius  form  of: 

kf  =  6.40  X  10'^  r-‘  exp(-38.370/r)  ■ 

mol  sec 

In  the  high  temperature  range,  the  exponent  on  the 
temperature  term  should  be  a  nonpositive  value.  It 
was  arbitrarily  set  to  -1  and  the  coefficient  was 
determined  via  a  curve  fit  to  existing  experimen¬ 
tal  data.  The  second  rate  is  one  determined  by 
Monat.^“  which  is  of  the  form: 
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kf  -  1.84  X  10'^  exp(-38,  370/T)  — 

mol  sec 

This  second  formulation  is  the  result  of  shock-tube 
studies  that  covered  a  temperature  range  of  2,000  K 
to  4.000  K.  The  temperatures  of  the  weak  bow-shock 
in  the  flow  are  of  the  order  of  10.000  K;  however,  the 
Monat  rate  is  at  least  derived  directly  from  experi¬ 
ment.  The  different  temperature  dependence  of  the 
two  rates,  especially  at  high  temperatures,  is  shown 
in  Figure  2.  For  these  reasons,  the  Monat  rate  is  a 
good  second  choice  to  determine  the  sensitivity  of 
the  predictions  to  the  exchange  reaction  rate. 

Radiance  Calculations 


NO  H-  O  ^  O2  +  N 
N2  -f  O  ^  NO  -h  N 

With  the  first  of  these  two  reactions,  the  back¬ 
wards  reaction  is  controlled  by  the  concentration  of 
atomic  nitrogen.  At  high  altitude,  the  dissociation 


In  the  final  step  of  the  solution  procedure,  the  flow 
field  solution,  with  correct  concentrations  of  both 
the  major  and  trace  species,  is  analyzed  with  the 
NEQAIR  code.  This  code  uses  input  from  a  calcu¬ 
lated  flow  field  and  can  be  used  to  compute  emis¬ 
sion  predictions  at  particular  wavelength  for  a  given 
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Figure  2:  Comparison  of  the  Park  and  Monat  Re¬ 
action  Rate  Coefficients  for  the  N2  +  O  — ^  NO  -hN 
Exchange  Reaction 

species.  Corresponding  to  the  data  that  will  be 
taken  during  the  Skipper  flight  experiment.  230  nm 
±  5  nm  NO  emissions  are  integrated  along  the  axis 
and  131.0  nm  atomic  O  emissions  are  integrated 
along  a  line  off  the  axis. 

Results 

At  all  altitudes,  three  separate  simulations  are  per¬ 
formed.  The  first  uses  the  MSIS  generated  atmo¬ 
sphere  and  the  Park  reaction  rate  for  the  second  ex¬ 
change  reaction.  Given  that  the  MSIS  atmosphere 
is  likely  the  more  accurate  prediction  of  the  actual 
atmospheric  structure  that  will  be  present  during 
reentry  and  that  the  Park  reaction  rate  is  applicable 
over  the  entire  temperature  range,  this  first  calcu¬ 
lation  is  considered  the  baseline  solution  for  each 
altitude.  In  order  to  investigate  the  effects  of  chang¬ 
ing  the  reaction  rate,  the  second  calculation  made 
at  each  altitude  uses  the  same  atmosphere,  but  em¬ 
ploys  the  Monat  rate  for  the  second  exchange  reac¬ 
tion.  Calculations  of  this  type  will  be  referred  to  as 
Monat  results.  The  third  calculation  at  each  alti¬ 
tude  uses  the  Park  reaction  rate,  but  employs  the 
standard  atmosphere  as  the.free-stream  conditions. 
These  calculations  will  be  referred  to  as  standard 
atmosphere  results. 

Nitric  Oxide  Radiance  Predictions 

The  230  nm  ±  5  nm  ultraviolet  radiance  predictions 
for  NO  at  all  altitudes  are  given  in  Figure  3.  The 
radiance  predictions  are  most  sensitive  to  the  NO 
number  density  along  the  stagnation  streamline  at 


Figure  3:  Nitric  Oxide  230  nm  ±  5  nm  Radiance 
Emission  Predictions 

the  point  of  the  peak  translational  temperature. 

At  100  km.  there  is  very  little  difference  between 
the  radiance  predictions  for  the  baseline  result  and 
the  standard  atmosphere  result.  This  is  expected  as 
the  NO  number  density  profiles  for  the  two  cases  are 
very  similar  along  the  entire  stagnation  streamline. 
The  locations  of  the  peaks  in  translational  tempera¬ 
ture  for  both  cases  is  relatively  unimportant:  regard¬ 
less  of  where  the  individual  peaks  are  located,  there 
is  little  difference  between  the  NO  number  densi¬ 
ties  computed  at  those  points.  \n  order  of  magni¬ 
tude  difference  exists  between  the  baseline  result  and 
Monat  result.  This  difference  arises  directly  from 
the  NO  number  density  profiles  for  which  there  is 
a  large  difference  at  the  point  of  peak  translational 
temperature. 

At  120  km.  the  radiance  prediction  of  the  stan¬ 
dard  atmosphere  calculation  is  about  a  factor  of  four 
higher  than  the  ba.seline  radiance  prediction.  This 
follows  directly  from  the  NO  number  density  profiles. 
A  difference  is  still  notable  between  the  baseline  pre¬ 
dictions  and  the  Monat  predictions.  Again,  this  is 
expected  from  examination  of  the  NO  number  den¬ 
sity  profiles.  The  overall  magnitude  of  the  difference 
is  less  than  is  the  case  for  the  100  km  case.  This  can 
be  attributed  to  the  reduced  magnitude  of  the  differ¬ 
ence  between  the  NO  number  densities  as  compared 
to  the  100  km  results. 

Relations  between  the  predictions  obtained  from 
the  various  140  km  and  160  km  calculations  are  qual¬ 
itatively  similar,  with  the  actual  differences  between 
cases  being  much  less  pronounced  at  these  higher  al¬ 
titudes  than  is  the  case  at  the  lower  altitudes.  At 
140  km,  the  difference  between  the  baseline  predic¬ 
tion  and  Monat  prediction  is  relatively  small,  due 
to  the  reduced  effect  of  chemistry  in  the  extremely 
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Figure  4:  Atomic  Oxygen  131  nm  Radiance  Emis¬ 
sion  Predictions 

rarefied  flow  field.  At  160  km.  the  insensitivity 
to  changes  in  the  reaction  rate  is  even  more  pro¬ 
nounced.  This  insensitivity  at  both  altitudes  carries 
over  to  the  standard  atmosphere  results.  Despite  the 
relatively  large  difference  between  the  MSIS  atmo¬ 
sphere  and  the  standard  atmosphere  at  the  higher 
altitudes,  little  difference  exists  between  the  predic¬ 
tions  for  these  two  conditions. 

Atomic  Oxygen  Radiance  Predictions 

Atomic  oxygen  131.0  nm  vacuum  ultraviolet  radi¬ 
ance  predictions  are  shown  in  Figure  4.  At  each 
altitude  only  the  results  from  the  baseline  and  stan¬ 
dard  atmosphere  calculations  are  plotted.  Since  the 
overlay  technique  only  changes  the  concentrations  in 
the  trace  species,  the  concentration  of  0  cannot  he 
affected  by  an  overlay  calculation.  Since  the  only 
difference  between  the  procedure  used  to  generate 
Monat  results  from  the  procedure  used  to  generate 
baseline  results  is  different  input  parameters  to  the 
overlay  calculation,  there  can  be  no  difference  in  the 
concentration  of  O  between  the  two  cases.  Thus,  the 
atomic  oxygen  radiance  predictions  of  the  Monat  re¬ 
sults  must  be  identical  to  the  baseline  predictions. 

Overall  trends  in  the  atomic  O  predictions  at  the 
various  altitudes  are  consistent  with  previous  results. 
The  difference  at  100  km  is  relatively  small,  again 
due  to  the  similarity  of  the  two  atmospheres  at  this 
altitude.  There  is  a  notable  difference  at  120  km. 
where  the  two  atmospheres  start  to  deviate.  At  the 
higher  altitudes,  the  difference  between  predictions 
is  not  cis  large  as  the  difference  between  the  two  at¬ 
mospheres  in  this  range  would  suggest.  Use  of  the 
different  atmospheres  results  in  differences  in  the 
temperature  and  atomic  O  density  profiles  along  the 


Figure  5:  XOo  732  nm  Radiance  Emission  Estimates 

line  of  integration.  These  differences  are  such  that 
they  counteract  each  other  to  some  extent,  resulting 
in  the  reduced  sensitivitv  of  the  predictions  at  the 
higher  altitudes. 

Visible  Emissions  Observed 
During  Orbit 

In  this  section  we  provide  estimates  of  visible  glow 
that  has  been  observed  previoush'  on  the  Atmo¬ 
spheric  Explorer  (AE)  satellite  and  on  the  Space 
Shuttle.  Orange  in  color,  this  glow  has  been  at¬ 
tributed  to  an  excited  state  of  XOo.  The  glow  model 
as  it  stands  has  too  many  undetermined  parameters 
to  allow  a  high-confidence  pre-flight  prediction  for 
the  Skipper  mission.  Therefore  we  have  estimated 
Skipper  intensities  using  an  extrapolation  of  the  data 
measured  by  the  Atmospheric  Explorer  satellite. 

The  extrapolation  is  shown  in  Fig.  5.  It  has  been 
assumed  that  the  Skipper  surface  temperature  is 
similar  to  that  of  AE  at  293  K,  The  dimensions 
of  the  AE  and  Skipper  spacecraft  are  also  similar. 
The  data  is  for  ram  measurement;  .\E  measurements 
were  at  an  angle  of  lOX  A  cosine  scaling  law  is  rec¬ 
ommended. 

The  peak  glow  intensity  is  anticipated  to  occur 
around  the  wavelength  of  732  nm  for  which  the  data 
is  presented.  Estimated  scaling  factors  for  other 
wavelengths  are  0.35,  0.6,  and  0.70  for  550  nm. 
600  nm.  and  800  nm  respectively. 

Xote  that  the  extrapolation  to  130  km  is  highly 
uncertain.  Skipper  measurements  in  this  lower  alti¬ 
tude  range  will  provide  a  major  contribuiou  to  our 
understanding  of  visible  glow  phenomena. 
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Visible  Shuttle  Glow:  Gleanings 


Presented  at  ARO  TI  Meeting 

BMDO  Missile  Signature  and  Aerothermochemistry  Program 

22  April  1996 


Visible  Shuttle  Glow:  Gleanings 
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The  Visible  Shuttle  Glow 


material  and  temperature  dependent 


Glow  Altitute  Dependence 


Atmosphere  Explorer  Glow  Data 
(Yee  and  Abreu,  1983) 


Brightness  (Rayleigh) 


EISG  Visible  Glow  Observations 
(Swenson  et  al.,  1995) 
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Transition  Fiow  on  Space  Shuttie-Diffuse  Refiection 

(BiRD,  1990) 


T ransistional  regime  effects  persist  to  about  300  KM 
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Collaboration  with  Prof.  Ian  Boyd,  Cornell  U. 


DSMC  Calculations 
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Vehicle  dimensions  (preliminary  caicuiation  for  AE) 
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CONCEPTUAL  MICRO-VIEW  OF 
ADSORBANT  BEHAVIOUR  RESULTING 
IN  VISIBLE  SHUTTLE  GLOW 
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Reaction  Cured  Glass  Coating 


C-8753 


Glow  Kinetic  Modeling  is  Complex 

96-1123 

e.g.,  surface  kinetics  for  O  +  NO 
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Thermal 

Desorption 
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Thermal 
desorption  rate 

Collisional 
desorption 
or  reaction 

0(g)  +  0(S) - >20(9) 

- >02(g) 

X(  g)  +  0(S) - >X(g)  +  0(g) 

0(g)  +  N0(S) - >0(g)  +  N0(g) 

- >N02(g) 

0(g)+N02(S)^0(g)+N02(g) 

N0(g)  +  0(S)^N0(g)+0(g) 
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X(g)  +  NO(S) - >X(g)  +  NO(g) 

Collisional 
desorption 
or  reaction 
efficiency 

Surface 

Reaction 

0(S)  +  N0(S) - ^N02(g), 

N02(S) 

0(S)  +  0(S) - >02(g),02(S) 

Reaction  rate 

*g  ^  Gas  Phase  Species 
'S  =>  Surace  Adsorbed  Species 


Now  add  other  mechanisms,  N,  NO,,  etc 
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Skipper  Glow  Predictions 
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Glow  Temperature  Dependence 
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EISG  data  shows  lower  AE 


PSI  Glow  Data  vs  Temperature  (Oakes) 
(800  nm,  aluminum  substrate) 
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